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o^ 

rv^ Abstract: We present a hierarchy of quantum many-body states among which 

many examples of topological order can be identified by construction. We de- 
'~~ ' fine these states in terms of a general, basis-independent framework of tensor 

O networks based on the algebraic setting of finite-dimensional Hopf C*-algebras. 

At the top of the hierarchy we identify ground states of new topological lat- 
tice models extending Kitaev's quantum double models [55]. For these states we 
exhibit the mechanism responsible for their non-zero topological entanglement 
5^ entropy by constructing a renormalization group flow. Furthermore it is shown 

Ch that those states of the hierarchy associated with Kitaev's original quantum dou- 

I ble models are related to each other by the condensation of topological charges. 

■"O We conjecture that charge condensation is the physical mechanism underlying 

Ch the hierarchy in general. 

1. Introduction 

^ 1.1. Background. Quantum states of many-body systems, for instance lattice 

C^ models of quantum spins, furnish a rich physical arena where simplified mod- 

OO els of condensed matter phenomena such as high- Tc superconductivity can be 

^N studied [2]. Such quantum systems are, however, difficult to solve numerically, 

^1 since the dimension of the corresponding Hilbert space grows exponentially in 

t^^ the number of components. 

^^ Two promising developments tackling this difficulty have occurred in recent 

^^ times. Firstly, the advent of experimental settings where lattice models can be 

". . implemented whose parameters can be controlled very accurately [S]. Secondly 

^ - and relevant for the present article - the understanding of the entanglement 
structure of models with local Hamiltonians. It is believed that ground states of 
such models generically obey an area law for the entanglement entropy |47[|17j . 

rrt This has justified the development of a number of approaches where the physical 
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states of the model can be described by so-called tensor networks whose set of 
parameters grows only polynomially in the number of components and which are 
therefore suitable for numerical simulations. Among these are the tensor product 
states (TPS) [SS], also known as projected entangled-pair states (PEPS) [i5] . 
which generalize the structure of matrix product states (MPS) [19j to higher 
dimensions, and the multi-scale entanglement renormalisation ansatz (MERA) 
based on renormalisation group ideas p51H5] . 

In 2003 Kitaev proposed topological quantum computation [26] : this is an 
approach to quantum computation which is based on quantum many-body sys- 
tems exhibiting topological order, i.e. systems that are effectively described by 
a topological quantum field theory [35] . Topological phases appear in condensed 
matter physics, for instance in the fractional quantum Hall effect [IS]; they can 
also be studied in the context of quantum field theories [3j . In Kitaev's proposal 
information is encoded in nontrivial loops of the underlying surface [15] and is 
processed by the creation, braiding and annihilation of topological charges j211 
[55]. The simplest models proposed for this purpose are Kitaev's original quan- 
tum double models, which are based on discrete gauge theories [3S], and Levin 
and Wen's string-net models pP. The latter are believed to describe all parity 
and time-reversal invariant topological phases on the lattice and are constructed 
from the idea that topological phases correspond to fixed points of renormalisa- 
tion group procedures. In Kitaev's models, the space of superselection sectors is 
the set of irreducible representations of the quasitriangular Hopf algebra D(CG), 
Drinfeld's quantum double [T6] of the group algebra CG describing the gauge 
symmetry. We will therefore call these models D(CG)-modcls. The simplest such 
model is the well-known toric code which is based on G = Z2. 

Recently, tensor networks for topological lattice models have been introduced 
and it is this development which forms the starting point for the current work. 
In [33] a PEPS ansatz is devised for a distinguished ground state of the toric 
code. Ground states of general quantum double models based on a group alge- 
bra can be described as MERA states [T], as can ground states of string- net 
models [30]. Other tensor networks for string- nets are given in [Tl]l23j. 



1.2. Aims. In the following we will outline the two main goals of this article: 
the extension of Kitaev's quantum double models and the derivation of tensor 
network representations for these generalized models. 

The first aim of this article is to extend Kitaev's quantum double models [5S] 
from the case of the group algebra of a finite group G to a finite-dimensional Hopf 
algebra with certain properties, as anticipated by Kitaev in 1263. This sheds more 
light on the Hopf algebra point of view Kitaev outlined originally, in particular, it 
clarifies how the structure maps of the Hopf algebra enter in the definition of the 
model and how the quantum double of the Hopf algebra arises in the description 
of the model's excitations. It turns out that an involutory Hopf algebra (i.e. a 
Hopf algebra whose antipode map squares to the identity) with an additional *- 
algebra structure is sufficient in order to define the generalized quantum double 
model properly. As we will show both requirements are satisfied if one chooses 
a finite-dimensional Hopf G*-algebra (also known as a finite-dimensional Kac 
algebra). In order to illustrate the construction of our generalized model we 
give an example of a finite-dimensional Hopf G* -algebra which is nontrivial in 
the sense that it is neither a group algebra nor the dual of a group algebra. 
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This extends well-known instances of the family like the toric code for G ~ 
Z2 or the model based on G = 6*3 which is universal for topological quantum 
computation [55] , 

The second aim of the present article is to derive tensor network represen- 
tations for the ground states of these models. As such representations allow to 
map out the phase diagram of the model and determine phase transitions be- 
tween topologically and conventionally ordered phases it is crucial to understand 
the properties of tensor networks describing states of the model deep within the 
different phases. In order to derive these tensor network representations for quan- 
tum double model ground states we introduce a novel diagrammatic technique. 
It turns out that the tensor networks can be formulated basis-independently and 
that furthermore a single distinguished element of the Hopf C*-algebra, namely 
its unique Haar integral, plays the key role in the construction. Other than that 
only the structure maps of the Hopf C*-algebra are employed in the definition 
of the tensor network states. 

In addition to the stated aims above we are able to extend these ground 
states to a hierarchy of tensor network states for each finite-dimensional Hopf 
C* -algebra. This hierarchy of states is characterized by different values of the 
topological entanglement entropy 7 as defined in 27,32J, and hence these states 
represent different instances of topological order, in other words different unitary 
modular tensor categories (UMTCs). The way this hierarchy depends on the 
Hopf subalgebras of the original Hopf C*-algebra points towards condensation 
of topological charges [4,5J. Furthermore we show how different isomorphism 
classes of finite-dimensional Hopf C* -algebras exhibit different mechanisms for 
the non-zero 7. In particular, we explain how the boundary configurations of 
a region differ between models based on a group algebra, the dual of a group 
algebra and non-trivial finite-dimensional Hopf C*-algebras. 



1.3. Y>{<CG) -models and their tensor network representations. We will now re- 
view the basics of Kitaev's quantum double model and its tensor network rep- 
resentations in a language appropriate for our later discussion. 

The quantum double model presented in [26] is a quantum spin model defined 
on any finite oriented graph F = (y, E, F) with vertices V , edges E and faces F 
which can be embedded in a two-dimensional surface. Throughout the article 
we will denote the set of all edges connected to a vertex s € V hy E{s) and 
similarly E{p) will stand for the set of edges that form the boundary of the face 
p E F. In order to construct the Hilbert space of the quantum double model 
one assigns to each oriented edge the local Hilbert space CG whose canonical 
basis is {\g) \ g G G}. Reversing the orientation of an edge corresponds to the 
map \g) 1— )• \g~^). The Hamiltonian is assembled from two sets of operators that 
act on vertices or faces of F respectively. At each vertex s G V one may orient 
the adjacent edges such that they point inwards. Then the vertex operator A{s), 
which acts on -E(s), is given by 

' ' geG 

where L{g) : CG — ^ CG, \h) 1— >■ \gh) denotes the left regular representation of G. 
Hence the action of A(s) is a simultaneous left multiplication at each incident 
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edge averaged over the group G. In fact, it projects onto the trivial representa- 
tion. Using the orientation reversal isomorphism one may orient the edges around 
a face p G F in counterclockwise fashion. One then defines the face operator B{p) 
by its action on basis states \gi, . . . , gr): 



B[p) 9r 
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Again, this is a projector which projects onto configurations with a trivial prod- 
uct of group elements around the face. Finally, the Hamiltonian reads: 



U 



sQV 



A{s) 



peF 



B{p). 



(3) 



This is a sum of mutually commuting projectors and therefore exactly solvable. 



In order to find a tensor network representation for particular states occurring in 
the quantum double model one can start from its inherent (group) symmetries. 
We can assign a representation pp of G to each face p £ F and associate the 
tensor 



13 Pp 
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with each oriented edge. Here Pp{g)ai3 denotes a matrix element of the repre- 
sentation Pp and the red dot represents the physical index g with its orientation 
inherited from the underlying graph edge. The Greek letters attached to black 
arrows are called virtual indices. The tensor for a reversed graph edge naturally 
reads 

P Pp a P Pp a 



9 :— 
\ f 

& 7 



(5) 



Reversing virtual arrows is defined by 

(i Pp a 



Pp 



Pp^ 



f 9 
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(6) 



where the dual representation p* of a representation p is given by the matrix 
equation p*{g) :— p^{g~^). Note that this is equivalent to a reflection of the 
tensor about its vertical axis: 



/3 Pp 



Pp P 
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^ / 

5 S 



g . 



(7) 



The tensor network obtained by contracting all virtual indices around each face 
represents what we will call a group tensor network state in the following. 
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For example, on a square lattice the tensor representing a vertex with a 
particular edge orientation is given by 



Q3 P2 



Pi P3 
OLi 



I I 



Q2 

Pi l3 



= YlPji9j ^ffj + l)fta, + i 



^P4 



(8) 



where the index j labelling graph edges is assumed to be cyclic. Quite similarly, 
the tensor corresponding to a particular face reads 



P3 
93 



P4 94 

/?4' 



P5 



92 P2 






(9) 



where Xp is the character of the representation p. 

In fact, both this vertex and face tensor are partially contracted tensor net- 
works and represent linear maps from the uncontracted virtual indices to the 
physical ones. So a piece of a tensor network with open indices can be regarded 
as a PEPS projection map P^I36j . If Vp is the vector space associated with the 
representation pp then this linear map P : Vp^Vp^Vp' ^Vpi — ?► CG is given by 



,<5eG 



(10) 



9.a> 



for a single edge tensor A{pp, pp>). 

On the contrary, a fully contracted tensor network is a complex number which 
equals the inner product between a particular basis state and the group tensor 
network state. If viewed as a linear map from basis configurations to amplitudes 
we will call it a tensor trace. 

Clearly, a group tensor network state only depends on the isomorphism class 
of the representations pp assigned to the faces. In general, there is vastly more 
gauge freedom for a tensor network to represent the same physical state, however, 
this does not necessarily respect the G-action. 



If we choose all representations Pp to coincide with the left regular one the 
resulting group tensor network state will be a ground state of the D(CG)-model. 
Since its Hamiltonian is frustration free (i.e. a sum of mutually commuting terms) 
this can easily be shown by considering independently the action of A{s) and 
B{p) on the partially contracted tensor networks in (Is]) and Q. 

As an example, a ground state of the toric code is given by the group tensor 
network for the group Z2 = {e, a} with all faces carrying the regular repre- 
sentation. In that case the PEPS projection map reads P = |e) {(p'^l {(f>'^\ + 
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\a) (■(/'+ 1 (V'"*"! with the Bell states 

!</>+) = J2 ^(e)"/3 I"' /^) = Ie, e) + \a, a) , (11) 

a,/3eZ2 

|^+) = ^ L(a)„^ |a, /3) = |e, a) + \a, e) . (12) 

We would like to remark that this projection map coincides (up to a trivial 
isomorphism) with the one given in 44 . 

Later in the article we will turn to the physical significance of group tensor 
network states arising from representations of G other than the regular one. 

In order to build some more intuition for the main part of this article we briefly 
review the structure of a local tensor in a group tensor network state. It is easily 
seen that the tensor in Q factorizes with respect to partitioning the virtual 
indices into the sets {a, /3} and {7,(5}. This implies that the PEPS projection 



map ( 10 ) can be rewritten as: 



P=Y.\9){M9)\{^P'i9-')\ (13) 

geG 

where a state associated to the virtual indices in face p is given by: 

IMa)) = J2 Ppi9)o.f!\a,P}. (14) 

a,l3eG 

This means that one can carry out the contraction of the tensor network in two 
steps. First one contracts each virtual loop separately and then one needs to glue 
these pieces together. In a group tensor network this glueing process is simply 
given by 

\9) ® \9'') ^ \9) • (15) 

It turns out that this step becomes rather nontrivial once we generalize the quan- 
tum double model to finite-dimensional Hopf C*-algebras. In fact, we anticipate 
that this glueing step is one of the crucial ingredients for deriving the hierarchy 
of tensor network states in this article. 



1.4- Outline. The article is structured as follows. In Section [21 we will provide 
a guide to the language of Hopf algebras, introducing the necessary notation 
needed in order to extend the quantum double models based on finite groups to 
Hopf algebras. Characteristic examples of Hopf algrebras relevant for this work 
are supplied, too. In Section l3] we then generalize Kitaev's quantum double con- 
struction from group algebras to Hopf algebras, thereby introducing the physical 
model which stands in the centre of this work. In Section |4] we will solve this 
model by providing a tensor network representation of its ground state. This 
representation only involves the canonical structures associated to the underly- 
ing Hopf algebra and leads us to propose a novel hierarchy of tensor network 
states based on Hopf subalgebras. In Section [5] we calculate both the entangle- 
ment entropy and the topological entanglement entropy for distinguished states 
in the hierarchy. In section [6] we provide concluding remarks and an outlook on 
future work. 
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2. The language of Hopf algebras 

2.1. Hopf algebras. If one sets out to find an algebraic structure for symmetries 
of linear spaces with tensor product structure, such as encountered in many-body 
quantum physics, one is naturally led to the concept of a Hopf algebra. This is 
because the axioms of a Hopf algebra directly allow for duals and tensor products 
of representations. As general reference for this section we recommend |25j . 

First and foremost a Hopf algebra H is a. vector space (over some field k 
which we will take to be C) equipped with some additional structure. On the 
one hand there is an associative multiplication fi of vectors which is linear in 
each argument, hence one may define it as the linear map n : H <S) H ^ H with 

/i(a; (^ y) = xy (16) 

such that 

{xy)z = x{yz). (17) 

This multiplication is accompanied by a unit rj which can be defined as the linear 
injection of scalars into H, hence rj : k -^ H where 

VW = Mh- (18) 

The element Ij/ € i/ is a left and right unit for the multiplication fi. At this stage 
H has the structure of an algebra. The multiplication encodes the composition 
of symmetry transformations in a representation space. 

On the other hand there is a dual notion to multiplication. This dual linear 
map is the comultiplication A : H ^f H ® H and is usually written in Sweedler 
notatiorQ as 





A{x)-- 


(x) 


(g>x"^'^x^'^ 

(x) 


®x(2). 








(19) 


This serves as a 
be coassociative: 


shorthand for sums of the form 


E, 


xf^<^ 


^'- 


It 


is 


required to 



(Z\ (8) id) o Z\ = (id ®A) o A. (20) 

Elements which satisfy A{x) = x ® x are called grouplike and collected in the 
set G{H). This comultiplication has a companion called the counit e : H ^ k 
which is required to satisfy the axiom: 

^e(x')x" = ^a;'e(x")=x. (21) 

(x) (x) 

In other words, it neutralizes the comultiplication. With this structure alone, 
the vector space H is called a coalgebra. The comultiplication encodes how 
symmetry transformations act on tensor products of representation spaces. The 
counit provides the appropriate notion of a trivial representation, or invariance 
under symmetry transformations. 



^ If there are more than three tensor factors in iterated coproducts such as A^^>{x) 
J^fa;) ^ ® ^ ® ^ ® ^ ™^ ™'" consistently use superscripts instead of primes. 
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If both the algebra and coalgebra structure are compatible then H is called 
a bialgebra. The compatibility axioms read: 

A{xy)^Aix)Aiy), (22) 

A{1h) = Ih®Ih, (23) 

e{xy)=eix)e{y), (24) 

e(lff) = Ifc. (25) 

Finally, if there is a linear map S : H ^f H satisfying the axiom 

^ a;'5(:E") = e(a:) Iff = ^ 5(a;') x" (26) 

(a;) (x) 

then H = {H; jj,,!]; A,e; S) is called a Hopf algebra. The map S is called the 
antipode and has the following properties: 

S{xy)^S{y)S{x), (27) 

S(1h) - Iff, (28) 

Y, 5(:E)'(»5(2;)"-^^(a;")®5(x'), (29) 

(Six)) (x) 

e{S{x)) - e{x). (30) 

In any finite-dimensional Hopf algebra the order of the antipode S is finite |40j , 
hence S is an invertible map. We will always make this assumption in the sequel. 
The antipode is used to define dual (or conjugate) representations. 

We may also define opposite multiplication /x°p and comultiplication Z\'=°p by 

IJ°'^{x(S)y) = yx, (31) 

Z\™P(x) =^2;"«)a;' (32) 

(x) 



relative to the multiplication (16) and comultiplication (19). Then the sets 



Z{H) = {xeH\yy£H: /i°P(a; (g) y) = fi{x (g) y)} (33) 

Cocom(iJ) = {xeH\ A''°P{x) = A{x)} (34) 

are called the centre of H and the cocommutative elements of H respectively. 
If Z{H) = H then iJ is a commutative Hopf algebra, in the case Cocom(i/) = 
i/ it is a cocommutative Hopf algebra. Furthermore it turns out that i?°P = 
{H;^°P,i];A,€;S^^) is again a Hopf algebra, called the opposite Hopf algebra 
oi H. 

One can now show that for a given Hopf algebra H — {H; /i, rj; A, e; S) with 
underlying vector space H the dual vector space H* has again the structure of 
a Hopf algebra. More precisely, 

H* = {H*;A^,e^;f,^,r^^;S^) (35) 
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with the structure mapf[j as indicated is called the dual Hopf algebra of H. By 
the same token, the opposite Hopf algebra H°p has the natural dual 

X = {H°Py^{H*;A'^,e^;{fi°Pf,r^'^;{S-Y)- (36) 

Finally, there is the notion of a (two-sided) integral yl in a Hopf algebra H. 
Such an element is defined via 

xA = Ax = e{x)A (37) 

for all X £ H. In view of the role of the counit as a trivial representation, 
integrals are thus invariant elements under multiplication. The dual definition 
of an integral F E H* can be phrased as 

^x'r{x") = ^r{x')x" ^r{x)iH (38) 

[x] (x) 



for all X E H. This turns out to be equivalent to (37) once applied to H* . 

Particularly important is the notion of a Haar integral. This is the normalized 
version of an integral which is guaranteed to exist uniquely for finite-dimensional 
Hopf C*-algebras. Those are the Hopf algebras we will be occupied with mostly 
in the discussion of the physical model. We will supply a precise definition of 
Haar integrals in Section [X2l 



2.2. Quantum doubles as bicrossed products. In a general Hopf algebra H noth- 
ing can be said about how A and A'^"^ are related to each other. In other words, 
H can be non-cocommutative in the most unpleasant way. Therefore we seek 
Hopf algebras which are almost cocommutativq^ in a certain sense: they are 
equipped with a so-called quasitriangular structure which controls the extent to 
which the comultiplication fails to be cocommutative. In [TB] Drinfeld gave his 
celebrated quantum double construction which produces such a quasitriangular 
Hopf algebra from a given Hopf algebra. 

At the heart of his construction lies the idea of a bicrossed product. For 
example, in the case of the semidirect product of two groups one has an action 
of one group on the other and this action is used to define the multiplication 
in the product group. Similarly, one may define a bicrossed product of groups 
where there is an additional backaction of the second group on the first one. It is 
this concept that can be generalized to Hopf algebras and will yield the quantum 
double. 

So given a Hopf algebra H with invertible antipode S (in particular, any 
finite-dimensional Hopf algebra) one can construct a new Hopf algebra D(iJ) = 
X t^ H, the quantum double of H, as the bicrossed product of X = (iJ°P)* and 
H. This means that as a vector space D(iJ) simply equals H* (E) H. In order 

^ For a linear map F : U ^ V the transpose map F"^ : V — > U* is defined as usual by 
{F'^{a),x) := {a,F(x)) for all a e V* and a: e U. 
•^ In physical applications this is related to the exchange statistics of particles in two dimen- 
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to define its multiplication we need an action oi H on X and another one vice 
versa. This bicrossed structure is given by the actions 

>:HxX^X, a[>/ = ^/(5'-i(a")?a'), (39) 

(a) 

<:HxX^H, a<f ^^f{S-^{a"')a')a", (40) 

(a) 

where the question mark denotes the argument of the function a\> f, i.e. in this 
particular instance 

ia>f){x)=Y,fiS'Ha")xa') (41) 

(a) 

holds for any x Cz H. We will use this notation frequently later on. 
One can then show that 

(/ <E> a){g ® 6) = ^ /^(^-^(a'") ?a') ® a%, (42) 

(a) 

together with the canonical comultiplication on the tensor product of Hopf alge- 
bras, defines a valid Hopf algebra structure on the vector space D(i7). Since both 
the embeddings ix ■ X ^ D(i/), f ^^ J ®Ih and in : H ^ D(iJ), a^^ lx®a 
are algebra morphisms this multiplication formula is actually already determined 
by the so-called straightening formulaj^ 

af := (Ix ®a){f® Ih) - Y. /(^"'(«"') ?«') «"■ (43) 

(a) 



It should be noted that both underlying actions (391 and (40) are derived 



from the adjoint representation which for a, x £ H is defined by 

ad(a)(a;) = ^a'a;S'(a"). (44) 

(a) 



2.3. Hopf *- algebras, li k = C then a Hopf algebra if may sometimes be equipped 
with a so-called ^-structure that will become most important for physical ap- 
plications. Namely, it will allow us to define Hilbert spaces and unitarity. For a 
general reference on Hopf ^-algebras see [55] . 

Firstly, a conjugate-linear map * : H ^ H which satisfies 

{x*y ^ X (45) 

(xy)* = y*x* (46) 

is called an involution and turns H into a ^-algebra. It follows naturally that 

Ih = lif ■ (47) 



* This should be contrasted with the canonical (uncrossed) multiplication in the Hopf algebra 
X (i) H. In this case we simply have af = fa for any f & X and a (^ H since multiplication is 
defined componentwise. 
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If an involution * is compatible with comultiplication 

A{x*) ^ A{xy (48) 

then H is called a *-coalgebra. Here the involution oi H (E) H is defined by 
{x (E)y)* — X* (E)y* . In a *-coalgebra one always has 

e{x*)^7{xj (49) 

where the bar denotes complex conjugation. 

A bialgebra H with an involution for which it is both a *-algebra and a 
*-coalgebra is called a *-bialgebra. 

Interestingly, if a Hopf algebra H also has the structure of a *-bialgebra then 
the interplay between antipode and involution is already determined: 

S{S(x*y)=x. (50) 

Consequently, H is called a Hopf *-algebra. Note that its antipode is always 
invertible (even if H is not finite-dimensional) . 

Finally, one can show that the dual of a Hopf *-algebra H is again a Hopf 
*-algebra with the involution given by 



fix) = f {Six)*). (51) 

2.4- Trivial Hopf algebras and their quantum doubles. Here we briefly discuss 
some cases of how Hopf algebras can arise from groups. In particular we focus 
on how a group itself can be directly understood in the language of Hopf al- 
gebras. This will be the key for relating Kitaev's original construction to our 
generalization. 

2.4.1. Trivial Hopf algebras. A finite-dimensional Hopf algebra is said to be triv- 
ial if it is a group algebra or the dual of a group algebra for some finite group G. 
Let {g\ g € G} be a basis of the group algebra CG with multiplication 

^{g®h)^gh (52) 

and unit ?7(lc) — e. Its comultiplication and counit are given by 

A{g)^g®g, e(.g) = Ic (53) 

and extended by linearity. Similarly, the antipode map is defined by 

S{g)=9-\ (54) 

It is easy to see that every group algebra CG is cocommutative, but in general 
non-commutative. Furthermore, one may endow CG with the involution map 

g* = 9-' (55) 

which is extended by conjugate linearity. This turns CG into a Hopf G*-algebra. 
Its Haar integral reads 

^-^\T.9- (56) 

' ' see 
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The dual (CG)* of a group algebra coincides with the space C*^ of linear 
functions from the group to the field of complex numbers. For the dual basis {Sg \ 
g e G} multiplication and unit are defined by 

^g^h = ^g,h5g, Vi'^c) ^'^^g = '^ (57) 

where 5g /,, denotes the Kronecker delta. The comultiplication, counit and an- 
tipode maps are given by 

^'y^a) ^^^u® Sy, e{Sg) = Sgie) = Sg,e, S{6g) = Sg-1 (58) 

uv=g 

for each g E G. One may check that every function algebra C is commutative, 
but in general non-cocommutative. When contrasted with the corresponding 
group algebra, this is not too surprising: since CG and C"-^ are dual to each 
other, one obtains the multiplication of one algebra from the comultiplication of 
the other one and vice versa. Again, C'' is a Hopf G*-algebra via 

S*g = Sg (59) 

and conjugate-linear extensiorr] . Finally, 

(^ = Se (60) 

is the Haar integral of C*^. 

2.4-2. Quantum doubles of trivial Hopf algebras. As a vector space the quantum 
double of a group algebra is nothing but D(CG) = C*^ ® CG with basis {Sx ® g \ 
{x,g) g G X G}. The fact that the quantum double is a crossed product is 
reflected in the multiplication map 

{Sx<»9){Sy<S)h) =6x,gyg-iSx<S)gh (61) 

where the algebra acts on itself by conjugation. Actually, this is fully determined 
by the straightening formula 

gSx -.^ {l<»g){Sx®e)^Sgxg-i<»g. (62) 

(and the embeddings Sx '-^ Sx'® e and g t-^ 1® g being algebra morphisms) Here 
gSx is a frequently used shorthand notation and should be clear from context 
whenever it appears. The unit of D(CG) is 1 (g) e. Comultiplication and counit 
are given by 

A{Sx®g)^Y.'^^'"'^9)®{Su®9), (63) 

UV — X 

e{5x (g) g) ^ e{g) Sx{e) ^ Sx^e (64) 

^ If G is an Abelian group then 5* = S„-i defines another Hopf *-structure for the standard 
convention. In fact, for G non-Abelian the choice of * or * corresponds to a choice of either the 
standard convention or the convention (a eg) b)* = 6* eg) a* for the action on tensor products. 
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on the basis elements. Note that the comultiplication in D(CG) is derived from 
the comultipHcation in CG and (C^)™? = ((CG)°p)* rather than from C^ itself! 
This is because as an algebra the quantum double is the bicrossed product 

D(CG) = ((CG)°P)*cx]CG. 
The antipode map is 

S{S^(S)g) = 5g-i^~.ig(»g~'^. (65) 

2.4.3. Hopf subalgebras of trivial Hopf algebras. The Hopf subalgebras of a Hopf 
algebra will be important for the construction of our hierarchy of Hopf tensor 
network states. 

li H = CG, all the Hopf subalgebras A and B can be characterized easily [4l|5]. 
One finds that the Hopf subalgebras of CG are precisely the group algebras CK 
foiKcGa, subgroup. Furthermore the Hopf subalgebras of C^ are exactly 
given by C*^/^ with N <\G a, normal subgroup. It is also not difficult to find the 
Haar integrals of these Hopf subalgebras. In the first case one has 

^K^W\Y.9- (66) 

' ' geK 

For the second case note that C*^/^ can be identified with the functions in 
C'^ which are constant on cosets of N. Now consider the character xn of G 
that factors through TV to the regular character of G/N. It is easy to see that 
when multiplied with the above functions xn acts as a two-sided integral. Hence 
taking normalization into account, one may convince oneself that \G/N\ xn 
is precisely the Haar integral of C'^/^ . By abuse of notation we may write it as: 



- E '^- (67) 



2.5. An example of a nontrivial finite- dimensional Hopf C* -algebra. We are pri- 
marily interested in finite-dimensional nontrivial Hopf G*-algebras, i.e. finite- 
dimensional nontrivial Hopf algebras over C with a compatible involution. In 
order to find an example of such a Hopf algebra it is enough to investigate 
semisimple Hopf algebras in low dimensions and check whether these support 
a compatible ^-structure. It turns out that any semisimple Hopf algebra H of 
dimension p or pq with prime numbers p < q is trivial [18 . The same applies to 
dimensions p^ 1341 . Hence one may expect nontrivial semisimple Hopf algebras 
only for dimH > 8. 

2.5.1. The Hopf algebra Hg. In fact, we need to look no further than dimiJ — 8 
for an example. At this dimension there exists a unique nontrivial semisimple 
Hopf algebra (up to isomorphism) which is commonly denoted by Hg and is due 
to Kac and Paljutkin 2^. Interestingly, it can also be endowed with a compatible 
involution which turns Hg into a Hopf G* -algebra. 

Following |33| this algebra can be presented by generators x, y, z and the 
relations 

x^^y^^l, z^ = -{l + x + y-xy), (68) 

xy — yx, zx — yz, zy — xz. (69) 
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It is now easy to see that B = {1, x, y, xy, z, zx, zy, zxy} is a basis of H^. We will 
denote the corresponding dual basis by {Si \ i e B}. The coalgebra structure is 
determined by 

A{x)=x<E)x, A{y)^y®y, (70) 

A{z) = -{l®l + y®l + \®x-y® x){z ® z), (71) 

e(x) = e(y) = e(z) = 1, (72) 

and from this it is evident that Q{Hg) ~ {l,x, y, a;j/} ~ Z2 x Z2. The antipode 
reads 

S{x) = X, S{y) = y, S'(z) = z. (73) 

Finally, the involution is given by 

X* = X, y* — y, z* = z^ = -{z + zx + zy — zxy) = z . (74) 

This shows in particular that antipode and involution act differently on the 
generators of Hg. 

The Haar integrals h G Hg, and (f> e H^ are given by: 

h = --{l + x + y + xy + z + zx + zy + zxy), (75) 

8 

cj)^Si. (76) 

Remarkably they look like the Haar integrals of a group algebra and its dual! The 
dual Haar integral (j) and the complex involution automatically yield a Hermitian 
inner product on ifg as discussed in Section |3.2[ Explicitly, we have for all 
a,b e Hs 

{a,b)^Si{a*b). (77) 

It follows that B is an orthonormal basis of Hs ■ 

2.5.2. Hopf subalgebras of Hg and its dual. The Hopf subalgebras oi H = Hg 
read {!}, {1, x} ~ CZ2, {1, y} ~ CZ2, {1, xy} ~ CZ2, {1, x, y, xy} ~ C(Z2 x Z2) 
and Hg,. This is because Hopf subalgebras of a Hopf algebra H correspond bi- 
jectively to sets of irreducible representations of H* which close under fusion 
and dualisation [31[S]. Now iJg has four one-dimensional irreducible representa- 
tions Qi, Qx, Qy and Qxy derived from Q{Hs) and the two-dimensional repre- 
sentation given by 

f-\ / z -\- zx z zx \ {^Q\ 

2 \zy — zxy zy -f zxy J ' ^ ' 

From the fusion rules 

Qg X Qg, = Qgg, (79) 

Qg X Q2 = Q2 (80) 

Q2 X Q2 = Ql + Qx + Qy + Qxy (81) 

for all g,g' € Q{Hg) one concludes that precisely the sets {Qi}, {Qi^Qx}, 

{Ql,Qy), {Ql,Qxy), {Ql,Qx,Qy,Qxy} and {Ql,Qx,Qy,Qxy,Q2} cloSC Under 
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fusion (and dualization) and therefore are the ones that correspond bijectively 
to Hopf subalgebras of Hg, . 

The Haar integrals of the Hopf subalgebras read 



^{1} == 1, 



Hx) 



Hv) 



^{xy) 



T-{x,y) 



(1 + a:), 

(1 + y), 

(1 + xy), 

(l + x + y + xy), 

(1 + X + y + xy + z + zx + zy + zxy). 



(82) 
(83) 

(84) 

(85) 

(86) 

(87) 



Since iJ| ~ i/g as Hopf algebras the Hopf subalgebras B C H^ are given as 
the images of the Hopf subalgebras of i/g under the isomorphism. A particular 
such Hopf isomorphism can be obtained as follows. From the dual basis {5a\ of 
the basis B C H^ define another basis {fa} C H^ via 



/l 


^ 5l+5x+5y+ S^y +S^+ S^x + S^y + 5y_^y, 


(88) 


/. 


= 5i+ 5x + 5y + 5xy — Sz — 6zx — Szy — Szxy, 


(89) 


fy 


= 5i-Sx - 5y + 5^y + iSz - iSzx - iSzy + iS^xy, 


(90) 


Jxy 


~ 5i — Sx — Sy + 5xy — i5z + i5zx + iSzy — i^zxy, 


(91) 


fz 


^5i+5x- 5y- 5xy +5z+ 5zx - Szy - Szxy, 


(92) 


Jzx 


^ Si- Sx+Sy - Sxy +Sz- 5zx + Szy - Szxy, 


(93) 


Jzy 


= Si + Sx — Sy — Sxy — Sz — 6zx + Szy + Szxy, 


(94) 


Jzxy 


= Si-5x +Sy - Sxy -Sz+ 6zx - Szy + Szxy 


(95) 



Then the map Hg — > Hg given by a h- > /„ for all a G i3 and linear extension is a 
Hopf isomorphism. 

Consequently, the Haar integrals (J)b read: 

^{fi} ^ Si+ Sx + Sy + Sxy + Sz + Szx - 

^{f=c) ^ Sl+ Sx+ Sy+ Sxy, 

1 + i 

'^{fy) = '^1 + Sxy H Tf- {Sz + Szxy) + 

1 - i 

'?^(/x„> = Sl+ Sxy H — (Sz + Szxy) + 

't'U^.Jy) =Sl+Sxy, 
= <5l. 



\- Szy + Szxy ~ e, 


(96) 




(97) 


2 [Szx+Szy), 


(98) 


1 +^ /r r ^ 
2 [Szx+Szy), 


(99) 




(100) 




(101) 



3. Constructing quantum double models from Hopf algebras 

The goal of this section is to construct a two-dimensional quantum spin model 
whose microscopic degrees of freedom are given by a finite-dimensional Hopf 
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algebra if, such that its emerging degrees of freedom are characterized by the 
quantum double 'D{H). This illuminates Kitaev's insight in [551. The model 
obtained in this fashion dynamically implements the quantum double construc- 
tion. This means that its Hilbert space aquires additional structure, namely that 
the graph underlying the spin model can be interpreted as a D(iJ)-module, or 
representation of D{H). State vectors describing (elementary) quasiparticle ex- 
citations above the ground state then live in the irreducible representations of 
D(if ) the Hilbert space decomposes into. In fact, for every irreducible represen- 
tation of the quantum double there will be a corresponding type of quasiparticle 
excitation. These quasiparticle excitations can then naturally be braided via the 
quasitriangular structure of D(ii) and, in general, exhibit exchange statistics 
beyond bosons or fermions. In order to render their topological nature manifest 
we would also like these quasiparticle excitations to be agnostic to the details of 
the microscopic background. In particular, they should only feel the topological 
properties of the underlying surface rather than the precise shape of the em- 
bedded graph. In other words, the quantum spin model we construct should be 
insensitive to the particular discretization chosen. It goes without saying that 
any discretization will typically have a fine granularity in order to make the 
condensed matter character of the model evident. 

Since there is a natural action of the quantum double D(ii) on the Hopf alge- 
bra H itself one can nevertheless regard a minimal graph consisting of precisely 
one edge as a representation of D(iJ) by identifying the edge with H. However, 
it is clear that this system cannot contain all irreducible representations of the 
quantum double. While this restriction can already be overcome by considering 
two graph edges associated with H ® H one still needs to go to a macroscopic 
regime from either of these small graphs. It is precisely the process of spatially 
extending the action of D(ii) from one edge to many edges which will eventually 
yield the Hamiltonian and thus define the quantum double model completely. 

In the following we will first construct a representation of D(iJ) from a small 
graph with just two edges. Next we will endow any larger graph obtained from a 
surface cellulation with a local D(iJ)-module structure. Finally we will introduce 
the Hilbert space for our generalized quantum double model and assemble the 
Hamiltonian from particular operators in the local representations of D(iJ). 



3.1. Graph representations of quantum doubles. Given a finite-dimensional Hopf 
algebra H and an oriented graph F ~ (V, E, F) we define the microscopic Hilbert 
space as Mr = 0eeB ^- ^^'" S*^*^ ^^ ^'^ ^"^"^ local operators Aa and Bf acting 
on Mr which represent both parts H and X = {H°^)* of the quantum dou- 
ble D(iJ) separately via the embeddings ix and in. Additionally we want these 
two operators to commute in such a way that they implement the straightening 



formula (431 of D(_ff). Put differently, we require the operators Aa and Bf to 
establish the bicrossed multiplication of the quantum double by interacting non- 
trivially on the intersection of their supports. Loosely speaking, this amounts to 
assembling the adjoint representation from smaller building blocks. Note that 
on the minimal graph comprising a single edge the adjoint representation can 
be turned into an operator "as is" . 

We begin with the small graph shown in Figure [T] and its associated Hilbert 
space Nr = F[®H . In the spirit of ^26] we define the following module structures 
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Fig. 1. Small graph representing H(^H. It affords a D{H)-module structure at the site {s,p) 
via the operators Aa{s,p) and Bf(s,p). 



on H for all a,x £ H and f £ X: 

L\{x):=ax, (102) 

U'L{x):=xS{a), (103) 

Ti{x):=Y,{f,x")x\ (104) 

TL{x):^Y.^f,S-\x'))x". (105) 

(x) 

More precisely, the operators L± define actions of H on itself while T± define 
actions of X on H . Furthermore these actions are intimately related: one may 
start with the left multiplication L+ and then canonically derive all other actions. 
As a consequence one has for example 

L^ ^S-oL^oS--!, (106) 

tL = SoTIoS-^. (107) 

This means that if one fixes an arbitrary pattern of edge orientations one may 
relate any other pattern to the original one using these relations. Unfortunately 
this does not treat all possible orientations on equal footing. We will resolve this 
issue in Section [3721 Also note that 

[Ll,Lt]^0, (108) 

[Tl,Ti]^Q (109) 

for arbitrary a,b £ H and f,gGX. 

Next we attribute operators Aa and Bf to the sites {s,p) of the graph F. 
Following l^nj , a site is defined as a pair consisting of a vertex s and one of its 
adjacent faces p. For the particular site {s,p) shown in Figure ll] we define the 
operators in terms of the above actions: 

Aais,p):=Y,L+^L''^\ (110) 

(a) 

(/) 

It is worth emphasizing that this definition only holds for the edge orientation 
chosen in the figure. However, using the conventions set forth by Figurel2]we can 
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I 

Fig. 2. A graph edge representing the Hopf algebra H via the H-module structures L±. At 
the same time it represents the Hopf algebra X = (i/°P)* v ia th e X- modu le structures T±. 
These are related to each other by means of the antipode in (|106| and (|107|l . 



extend this definition to arbitrarily oriented edges. Also note an important detail: 
our seemingly similar notation refers to two rather different comultiplications, 
namely, the one in H and the other one in X. 

We are now prepared to analyze the commutation properties of these opera- 
tors which we take both at the same site {s,p) for the moment. It is clear that 
the support of both operators coincides on this small graph, hence we can study 
the interaction of Aa and Bf without any additional complications. The follow- 
ing lemma shows that indeed the two operators define a representation of D(iJ) 
on H(g) H. 



Lemma 1. For the operators Aa{s,p) andBf{s,p) the straightening formula (43) 
holds in the form: 

^aBf - Y. Bf(^s-Ha"')^a')Aa"- (112) 

(a) 

Proof. The statement is proven by evaluating both sides of the straightening 
formula on arbitrary elements a:i,X2 € H. To this end consider first: 

AaBfixi ® X2) ^AaY, TL' {xi) ® TI" {X2) 
if) 

= AaY.Y.^f\S-\x',)){f",S-\x'2))x'[®x'i 

{^^) if) 

(^.) (/) (a) 

= Y. {f.S-\x',)S-\x',))a'x'{®x'iS{a") 

{a){xi) 

= Y {f,S^^{x[x2))ax'l(g)X2S{a"). 

(a)ix,) 

Note that in the third line we have used the opposite multiplication p°^ which 
is the appropriate dual to swap the comultiplication 

^/'®/" = (M°Pf(/) 
(/) 

in X from / onto its argument in H. 
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In a second step compare this with: 

^Bf(^S-Ha"')?a')Aa"{xi (» X2) = ^ B j.-(^g-^^(i)y,^(i)-^{a^^^ Xi (g) X2 S {a^^^)) 
(a) (a) 

(a) 

where we used the temporary abbreviation / = f(S^^{a^'^') 7a^^'). Now we have 
for the right-hand side 

(a)(/) 

= T. T. {f',S-'{{a(''^x,y)){a^'h,)"^ 

(a)(/)(a(2)xi) 

Y: (/",5-i[(.2 5(a(3)))'])(x,5(a(3)))" 
-E E E {LS-'[ix2Sia^^^)y]S-\iai^^x,y)) 

(a) (a(2)a;i) (a;2S(a(3))) 

= E E E {lS-'[{a(^^x,y{x2S{a('^)y])ia('^x,r^{x2S{a(^'>)y' 

(a) (a(2)xi)(2;2S(a<3))) 

and upon restoring / this becomes 

E E E {LS-\a(^^)S-'[{a('^x,y{x2S{a(^^)y]a(^^) 

(a) (a(2)xi) (a;2S(a(3))) 

= E E E (/,^-M(«^^'-i)'(-2 5(a(^)))V4^]a«) 

(a) (a(2)xi)(x2S(a(3))) 

= E E E (/,^-^((«<^^M45(a(^^)'«('')«w) 

(a)(x.)(a(2))(S(a(3))) 

(a(2))'V/®4'5(a(3))" 
= E E (/,5-i(a(2)x'i45(aW)'a(5))a«)a(3):z;'/®x^'5(aW)" 

(a)(x.)(S(aW)) 

= E E(/'^"'[«^'M4^((«^'^)")«^'^]«^'')«''^^1«4'^((«''')') 
(a)(x,)(aW) 
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= J2 f{S-\x\x'^) e{a')) a"x'l (g> x'i S{a"') 

(a)(x.) 
(a)(x.) 

Note that in the third Hne we have employed the fact that A is an algebra 
morphism. Since the above is true for all xi,X2 € H we have just proven the 
claim. D 



In order to obtain local D(_ff)-representations at the sites of arbitrary graphs we 
need to extend the above actions of Aa and Bf to larger vertices and faces. It 
will inevitably happen that the two operators Aa and Bf act on different sets 
of edges which do not fully coincide. We need to ensure that a) both operators 
continue to represent their respective parts H and X individually and b) the 
commutation relation arising from common edges still implements the bicrossed 
structure of the quantum double D(_ff). 

The extension can be done as follows for general vertices s and faces p. 

Definition 1. Let {s,p) be a site of the graph F, a €z H and f G X. Then we 
define vertex and face operators via 



(a) 



'Lf\xr.) 



E 

(a) 



-i a(i'a;i , 



(113) 



Bf{s,p) -A o 



TL {xi)®---®Ti (Xr) 

if) 



J2f{S-\x[)---S-\x',))< 



■^2 (114) 



Note that for Aa{s,p) the vertex s denotes the centre of the loop on the dual 
graph and p denotes the starting point within the loop for comultiplication. 
Analogously, for Bf{s,p) the starting point of the loop around the face p is 
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a;i 


P 


s 


' Xi 







Fig. 3. Graph used in the proof of Theorem 111 The commutation relation and hence the 
structure of T)(H) is determined on the intersection of the supports of Aa(s,p) and Bf{s,p). 



marked by s. For different edge orientations the actions L+ and T_ in (|113|) and 



(114 1 may need to be replaced per edge by L_ and T+ according to Figure [2J 

The idea behind this is: first one fixes the action oi a E H such that the 
coproduct of a is apphed to the edges of the vertex s in counterclockwise order. 
This is indicated by the arrow winding around s. When defining the action of 
/ e X on a face it then appears at first there might be two choices for the 
orientation of the coproduct of /. It turns out that only one of them, namely 
distributing the coproduct of / in cZocfcwise orientation on the graph, will yield 
a meaningful theory. The orientation of the coproduct is again indicated by the 
arrow winding around p. 

At this point the reader may wonder why one does not take the function / 
directly from H* (the dual of H without flipped comultiplication) and the fol- 
lowing more symmetric definition of the action (114): 



Bf{s,p) -r 



rj 



if) 



>T!!^\xr) 






(115) 



Note that here the arrow winds around p in counterclockwise orientation since 
/ G H* has a different coproduct now as compared to (114). Clearly, both (114) 



and (115) define the same (algebra) action on the boundary edges of the face. 



However, as cumbersome as taking f £ X in the definition of Bf{s,p) may seem 
it will make the next task tremendously easier: making contact with the Drinfeld 
double T){H) — X \xi H which these operators Aa and Bf represent. 



Theorem 1. Let H a finite- dimensional Hopf algebra. Then each site {s,p) of 
the graph F supports a D{H)-rnodule structure given by the operators Aa{s,p) 
and Bf{s,p) from DefinitionUl 

Proof. In order to keep the notation simple we only consider sites such as the 
one depicted in Figure |3] The generality of the argument will not be affected. 

In order to prove the D(iJ)-module structure at site {s,p) it is again enough 
to show that the straightening formula holds. Let a G H, f G X and evaluate 
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AaBf on arbitrary edges Xi Cz H: 

AaBfixi (E> ■ ■ ■ (» xe) ^ AaY, fiS'H^i) •5"' (4) -5"' (4) ^"'K)) 

x'( (8) X2 (8) X3 (8) X4 (8 Xg 8) Xg 

=. J2 f{s-H^4) S-H4) s-'i^'e) s-'i^'i)) 

^ a(i)x'/ (8 X2 5(a(2)) a^^^xg (8 x^' S{a'^^^) (8 x^' (8 x^' 

(a) 

= ^(/,5-i(x;x^x^xl)) 
(a)(^.) 

a^i^x'/ (8 X2 S-Cflf^)) ® a(3)x3 (8 x^' S'(a('^)) (8 x^' (8 x^'. 



Compare this with the case: 

5Z ^/(S-i(a"')?a')^a"(xi (8 • • • (8 Xg) 
(a) 

" X! ■S/(S-i(a(6))?a(i)) a^^^xi (8 X2 S{a^^^'>) (8 a^^^^xg 8) X4 S'(a(^)) (8 xg (8 Xg. 

(a) 

Again, abbreviating / = /(5'^^(a*^^'') ?a(^^) we obtain for the right-hand side 
T. T. T. T. {.f,S-'[{x,S{a('^))']S-\x',)S-\x',)S-\{a('^x,y)) 

{a(-^^x,y (8 X2 5(a(3)) ® a^^^Xg (8 (x4 ^(a^""^)))" (8 x^ (8 x^' 

= EE E E(/>"'[45((«(^^)")]^''(4)5"^(4)^-^((«^'^)'4)) 

(o) (a;i) (a(2)) (a(5)) 

(a(2))"a;'/ ^ a;2 ^(a(3)) ^ ^^^^xg (8 x:;' ^((a^^))') <8 x'5' (8 x^' 

= E E E(/'^"M(«^'')'4444^((«^'^^)")]) 

(a)(2;.)(a(2)) (a(5)) 

(a(2)) V/ ® X2 5(a(3)) ® flWxg (8 x'l ^((a^^))') <8 x'g' (8 x^' 
and now restoring / this becomes 
^(/(5-i(a(«))?a(i)),5-i(a(2)x;x^x^x;5(a(^)))) 

(a)(a:,) 

a(3)x'/ (8 X2 S'Caf'^)) (8 a^^^xg (8 x^' S'(a(s)) (8 x^' (8 x^' 
= ^ (/,5-i(a(«))a(^)5-i(x;x^x^xl)5-i(a(2))aa)) 

(a)(x,) 

a'^^x'/ ® X2 5(a('')) ® a^^^xg ® x^^ S'(a(^)) ® x'5' x'g' 
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= Y,{f,S-\x\x',x',x',)) 

(a)(x.) 

a^^'^x'i (E) X2 5(a(2)) (g, a^^'fxs (g) 4 ^(a^^)) (g x'5' (g <. 
This proves the statement. D 



3.2. Hilbert space. In order to obtain a physical system from the above discussion 
we need to define both a Hilbert space and a Hamiltonian for our topological 
lattice model. In particular, we need to find Hopf *-algebras H that allow for an 
inner product and *-representations. This means that adjoint operators (which 
as usual are defined relative to the inner product) in a representation of H are 
compatible with the ^-structure of H itself. It turns out we can reach both 
goals by requiring i/ to be a finite-dimensional Hopf C* -algebra. Such a Hopf 
algebra comes endowed with a unique element called the Haar integral which 
will naturally define both the inner product and the Hamiltonian. 

In order to define the Hilbert space we begin with the following proposition 
which is obtained from [571I2U] . 

Proposition 1. Let H be a finite- dimensional Hopf C* -algebra. Then S"^ = id 
and there exists a unique two-sided integral h ^ H with the following properties: 

1. h^ = h, 

2. h* = h, 

3. S{h) = h, 

4. he Cocom(iJ). 

Furthermore, H* is a Hopf C* -algebra again and its unique integral (f> G H* 
satisfying I-4 is a faithful positive functional, or trace, on H . 

As a first consequence we can easily resolve the issue of edge orientation: since 
the antipode is now involutive we define the reversal of an edge e simply by 

Xe ^ S{Xe). (116) 



This is obviously compatible with the actions (102), (103), (104) and (105). 
At the same time we no longer need to pick a distinguished pattern of edge 
orientations, rather all patterns are equivalent. 

Although the preceding proposition tells us that is a positive trace on H, 
we need to know its precise relationship with the usual trace tr//(a) = tr(L'^) 
found in the literature on Hopf algebras. Setting \H\ := dimH we have 

Lemma 2. Let H a finite- dimensional Hopf C* -algebra. Then 

tTH = \H\-cj) (117) 

holds where (p G H* is the Haar functional on H. 
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Proof. By the above Proposition we know that S^ ~ id. In this case we have 

tTH = e{A)cj, (118) 

for the integral A ^ H which is normahzed such that 4>{A) = 1 [3D]. Since 
(f){h) — \H\~ for the Haar integral h ^ H [29j we actually have that A = \H\ h. 
Finally, e{h) — 1 concludes the proof. D 

Since is a faithful positive trace on H we can derive a Hermitian inner 
product on H from it by setting 

{a,b)H = <l>{a*b). (119) 

This inner product now turns both the module structures L± and T± on H 



into ^-representations. Indeed, relative to (119) the adjoint map of Lj. is given 
by {Lj.y — Lj. because of 

{x,Ll{y))^^{x*ay) = (0, (a*x)*6) = {Lfix),y). (120) 

An easy but tedious calculation shows that (T^)^ = TJ. holds, too. Remember 



that /* is given by (51 1. 

This means that for the operators Aa and Bf which represent the Drinfeld 
double D(iJ) the adjoint operators are given by: 

Al{s,p) = Aa.{s,p), (121) 

B}{s,p)^Bf,is,p). (122) 



3.3. Hamiltonian. It remains to specify a Hamiltonian for the model. In analogy 
to [2B] we would like to get a frustration-free Hamiltonian, i.e. a sum of com- 
nuiting terms, and we would like to derive it from the local operators Aa and Bf 
defined previously. Hence we need to identify a subset of these operators such 
that they mutually commute with each other. 

Before anything else it is natural to analyze the commutation relation between 
Aa and Bf at the same site {s,p) of the graph F. Suppose in the following that 
both a € H and f € X are such that {fg,a) — {gf,a) for any f,g ^ X and 
f{xy) = f{yx) for any x,y € i/PlThen at the level of D(i?) we have 

(a) 

= J2f{a'S-\a"')7)a" 

(a) 

= E/(«"'^"'(0?)«' 

(a) 

= J2f{eia")?)a' 

(a) 

= fa, 



^ It is easy to see that this is equivalent to a £ Cocom(H) and / £ Coconi(X). Clearly, 
Coconi(_ff) and Cocom(X) are defined relative to the respective coproducts in H and X, 
which are fundamentally different. 
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where we used the cocommutativity of / (or a respectively) in the second (third) 
hne and the skew-antipode S~^ in the fourth one. In other words, the straighten- 



ing formula (43) becomes trivial for such elements. Since for a fixed site (s,p) the 
operators Aa{s,p) and Bf{s,p) form a representation of D(_ff) this commutation 
relation immediately carries over to: 

Aais,p) Bfi.S,p) = Bfis,p) Aa{s,p). (123) 

Interestingly, with this choice of a and / the full specification of a site {s,p) 
becomes unnecessary for the operators Aa and Bf. Indeed, since coproducts can 
now be permuted cyclically, the starting point p of the dual loop used to define 
the vertex operator Aa{s,p) loses any significance. All that matters for the action 
is the vertex s itself. In order to reflect this fact we set for all a E Cocom(iJ): 

Aais):^Aa{s,p). (124) 

By the same token a face operator Bf{s,p) only cares about the face p once we 
choose / e Cocom(X): 

Bf{p):^Bf{s,p). (125) 

Having restricted the possible candidates for the terms in the Hamiltonian 
by exploiting the D(_ff )-niodule structure at a single fixed site we need to ensure 
subsequently that all vertex operators Aa{s) commute among themselves, too. 
Observe first that 

[Aa(s),Ab(s')]=0 (126) 

for any a,b E H whenever two vertices s and s' do not coincide. Indeed, if Aa{s) 
acts on a common edge via L+ then Af,(s') acts on the same edge via L_ and 



vice versa, hence both operators commute by (108). If there are no common 
edges then there is nothing to show. On the other hand if s = s' then the edge 
set E{s) has an _ff- module structure (inherited from the D(_ff)-module structure 
at any site that contains s) which implies 

[Aa{s),At{s)]^A[aMis). (127) 

This suggests to further narrow down our set of candidate vertex operators Aa{s) 
by additionally requiring a G Z(i/) where Z(iJ) denotes the centre of the Hopf 
algebra H. By analogy, we are naturally led to assume / £ Z(X) n Cocom(X) 
in the following if we want all face operators Bf(p) to commute. 

Finally we would like to remark that Aa{s) and Bf{p) trivially commute if 
the pair (s,p) is not a site. In fact, the sets of edges they act on are even disjoint 
in this case. 

Recalling the properties of the Haar integral from Proposition [T] we now state 
the main result of this section. 

Theorem 2 (Generalized quantum double model). Let H a finite- dimensional 
Hopf C* -algebra with Haar integral h and Haar functional (j) and let F a graph. 
Furthermore for each s eV and p E F define the projectors 

A{s) -.^ Ah{s,p), (128) 

Bip) := B4s,p). (129) 
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Then 



n = -Y,A{s)-Y,B{p) 



(130) 



is a local, frustration- free Hamiltonian defining the D(H)-model. 

Proof. First observe that 

Aisf = A^,2{s,p)^Ais) (131) 

by Proposition fl] The same argument shows that B{p) is a projector, too. 

By the preceding discussion and PropositionII|it is clear that all lo cal terms A{s) 
and B{p) commute with each other. Furthermore, ( |121[ ) and (122 1 imply that 
they are Hermitian. D 

For the reader's convenience we give a short summary of how these operators 
act on the graph: 



A{s) ...- 



B{p) 



-001 = ^2 •■- 

(h) 



h'-^h 



-hi^^x, , 



(132) 



/.<'■) a 






Alternatively, for a different orientation of the graph edges the face operator acts 
as follows: 



Bip) 



(^0 



[xi ■ ■ ■ Xj.) 




(134) 



It is no coincidence that already at this stage do the Haar integrals of H and 
H* reveal themselves as the crucial ingredients for the generalized quantum dou- 
ble models. Since we set out to construct a quantum spin model whose (elemen- 
tary) quasiparticle excitations are characterized by irreducible representations 
of D(iy) the ground state sector with no quasiparticles present necessarily has 
the structure of a trivial representation of D(iJ) locally. Precisely this is what 
the element (j) ® h € D(^) embodies. As the Hamiltonian intrinsically encodes 
information about the ground state sector it should not surprise the reader that 
the Haar integrals h and <j) play such a prominent role. 
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3.4- Comparison with D{CG)-models. We would now like to briefly comment 
on the relation between our generalized quantum double models and Kitaev's 
original construction in [261 . 



First note that in the case H = CG the terms given by ( 132 ) and (133) in 
the above Hamiltonian reduce to the operators of (II]) and ([2]) which are the 
ones Kitaev employed for his quantum double models based on a group G. This 
becomes clear from Section [2.4. 1| where the relevant expressions for h and </> are 
listed. 

Secondly, the entire theory of ribbon operators readily carries over from |26) 
to our generalized quantum double models. This is because ribbon operators are 
constructed from certain elementary operators associated with two types of tri- 
angles any given ribbon path decomposes into. These operators are nothing but 
the L± and T± that implement the H- and X-module structures. Furthermore 
patching together a ribbon operator from those elementary pieces only involves 
the structure maps of the Hopf algebra H itself. Actually, we even used rib- 
bon paths in order to define the local D(iJ)-module structures (and hence the 
Hamiltonian) without saying so. 

Finally, our generalized quantum double models inherit all the beautiful topo- 
logical properties of the original since these follow exclusively from the algebraic 
structure of ribbon operators. In particular, these features include the degener- 
acy of the ground state sector as well as the exotic statistics of the quasiparticle 
excitations whose anyonic nature is revealed via braiding and fusion operations. 



4. A hierarchy of topological tensor network states 

In this section we develop a general diagrammatic language for tensor network 
states built from finite-dimensional Hopf C*-algebras. Underlying surfaces both 
with and without boundaries are considered and we show how to naturally de- 
scribe subsystems. Using this framework we solve the generalized quantum dou- 
ble model introduced in the preceding section by providing a tensor network 
representation for one of its ground states. Any other energy eigenstate can be 
obtained from there by an appropriate ribbon operator. The tensor network rep- 
resentation for that ground state only involves the canonical structures of the 
underlying Hopf C*-algebra H and its dual: multiplication, comultiplication, 
antipode and Haar integral. This insight leads us to propose a novel hierarchy 
of tensor network states based on Hopf subalgebras. For H — CG we are able 
to completely classify this hierarchy of states in terms of charge condensation. 
Finally we describe the relation between our Hopf tensor network language and 
the usual formulation of PEPS. 

Unless otherwise noted, from now on H will be a finite-dimensional Hopf 
C*-algebra with Haar integral h £ H and Haar functional £ H* . 



4.I. Tensor traces. As we discussed in the introduction, the fully contracted 
tensor network (which is a complex number) for a certain ground state of the 
D(CG)-model on the graph F can be interpreted as a collection of virtual loops 
in the faces of F that are suitably glued together to form the physical degrees 
of freedom. 
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We can extend this idea to the case of any finite-dimensional Hopf C*- 
algebra H now. In each face p d F place a virtual loop oriented in counterclock- 
wise direction and associate a function fpGXto this loop. For the moment we 
may restrict to fp G Cocom(X). With each oriented edge e G -E we associate an 
algebra element Xe (z H which splits into two parts as follows: 

((5'(g)id)oZ\)(a;e) =^5(x^)®<. (135) 

Subsequently, we attribute a;" to the left adjacent face of e and S{x'^) to the right 
one. A virtual loop in face p is then evaluated by taking the clockwise product 
of all elements thus associated with the loop from the surrounding edges of p. 
The result is then fed to the function fp. 

In order to simplify arguments we introduce a diagrammatic notation which 
encodes calculations with these tensor networks. In diagrammatic language the 
evaluation rule just described reads 



Y,fp{Six'J...)f,{x':...). (136) 



Since fp is assumed cocommutative its argument can be permuted cyclically and 
we may start both virtual loops around p and q at the edge e without loss of 
generality. The clockwise order of the product remains important though. Both 
the dashed lines and the ellipses denote the remaining degrees of freedom of the 



faces p and q respectively. Note that the glueing procedure generalizing (151 is 



implemented by the coproduct. We will introduce a full description of this later 



in ( 143 ) and ( 144 ) . Finally, a reversed edge is resolved via 

Jp ' ' Jp 



^..:=^5(..). (137) 

1 /, ; ; /, ; 

In summary we have the following 

Definition 2 (Hopf tensor trace without boundary). For each p G F let 
fp G Cocom(X). Then the Hopf tensor trace associated with the graph F is the 
function ttrr : i/^l^l X^l^l ^ C, 

(g)xe(g)/p>->ttrr({xe};{/p}) (138) 

eeE peF 



which is defined via diagrams and the evaluation rules (136) and (137) 



From this definition it is clear that the tensor trace is linear in each argument. It 
can be regarded as the wavefunction amplitude of a quantum many-body state 
as we will see shortly. 
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Fig. 4. Ordinary faces vs. boundary "faces" . Wherever the boundary of the surface M does 
not coincide with an edge of the graph F it is displayed in grey. Upon smooth deformation of 
this boundary the grey vertices on either side can actually be identified with each other. This 
modifies F and thus creates new faces whose preimages on the left-hand side we call boundary 
"faces". Boundary edges are drawn in red. 



In the following we would like to investigate Hopf tensor networks on graphs 
which are embedded in a surface M with boundary as well as Hopf tensor net- 
works at interfaces between regions i? of a surface with or without boundary. 

In general the edge set E = Em of the graph F embedded in M may be 
decomposed into the disjoint union of interior edges Em° and boundary edges 
EdM (see Figure l4|. On the other hand, the set of faces F = Fm = Fm° only 
contains interior faces by definition. However, one may collect incomplete faces 
(where two vertices of degree 1 can be connected/merged to yield a genuine 
new face) into the set Fqm of boundary "faces" as illustrated in Figure |4] By 
construction these are not faces of the original graph, hence F n Fqm = CK 

Furthermore we may collect all boundary segments into the set Xqm = Eqm^ 
FdM- Once we fix a distinguished segment on the boundary this set is equipped 
with a natural ordering inherited from the orientation of dM . Without loss of 
generality we will take any boundary to be oriented in counterclockwise direction 
with respect to its interior for the remainder of the article. 

In order to take nontrivial boundaries into account we need to refine our 
diagrammatic notation for tensor networks. Pick an arbitrary edge e € E. Since 
it belongs to the boundary of at least one face p in either F or Fqm (i-e. p is 
either an interior or a boundary face) we may define the following elementary 
diagram for any elements /„ € X and x^ € H simply as their canonical pairing: 



fp 



^ ,. 



• Jp\^e)' 



(139) 



Different orientations of graph edges and virtual loops are resolved via 

fp jp 

^ ^/^ t^ S(.;e) 

fp s^iu) 



(140) 



where all these elementary diagrams are assumed to be invariant under arbitrary 
rotations, for instance: 

" = A (112) 



-A. 
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Note also that S'^ = {S^'^V is the correct antipode oi X ^ {H°p)* . As a 
consequence, any of the above elementary diagrams has the same value as its 
mirror image under reflection about a vertical axis. In fact, (1401 and (141) are 



mirror images of each other in that sense. Furthermore for some fp (such as the 
dual Haar integral (j)) it may happen that S'^{fp) — fp and the loop orientation 
may become unimportant. 

If the face p has edges other than e in its boundary we may extend the above 
diagrams as follows. Pick another edge e' G E{p) which shares a common vertex 
with e. Then for any Xe' G H we define a "virtual" glueing operation by 







Jp\^e'^e) 



(143) 



where the arrows indicate the order in which the coproduct of fp is applied to 
the elementary diagrams. The black dot denotes the origin for this comultipli- 
cation. While this origin is very important in general it can be neglected iff fp 
is cocomniutative. In that case we will simply omit the corresponding dot from 
the diagram as we did implicitly in Definition [2] In any case one needs to pay 
attention to the correct comultiplication in X which causes the product around 
the edges of p to be taken in clockwise order. 

Finally for any interior edge e G Em° with adjacent interior or boundary 
faces p,q G F U Fqm we pick Xe G H and fp.fq G X arbitrarily and define a 
"physical" glueing operation by 



(-) ^ '" 



fr. 



(144) 



where the order of comultiplication is determined by the orientation of the un- 
derlying graph edge. Consequently one has for instance 



fv 



fp 



^a. 



fp 








f I f > 

fp 

1^ m—, S(x,)" 
(Six,)) <^ 



s(x^) (145) 



which is perfectly compatible with (1161 as expected. In contrast to elementary 



diagrams such a composite diagram may only be invariant under rotation by vr 
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Fig. 5. Diagram encoding the tensor trace ttr({a;e}; {/p}; {Ve}, {Sq})- While the interior de- 
grees of freedom {xe} C H and {fp} C X are only shown partially the boundary degrees of 
freedom {j/e} C H and {gq} C X are labelled in such a way that the ordering of the boundary 
is evident. 



if S{xe) = Xe- Also note that in general a simultaneous reversal of both virtual 
loops is not given by reflecting a composite diagram about its vertical axis: 



jp jp 

(-.) <^ "^ (-e) <^ "" 



(146) 



Jp r ,■ Jp 

Rather, equality holds iff Xe G Cocom(i7). This should be compared with (It]). 



Thus by startin g fro m the elementary diagram ( 139 ) we have reexpressed the 
evaluation rule ( 136 ) for an interior edge entirely in terms of "virtual" ( 143 1 and 



physical" (144) glueing operations. These are given by comultiplication in X 



and H respectively. This means that an interior edge is formed by appropriately 
glueing two boundary edges together via comultiplication. Alternatively, one 
may regard this as glueing together two virtual interior loops. Virtual loops 
themselves are assembled via glueing together smaller loop pieces. 
In summary we have the following general 

Definition 3 (Hopf tensor trace). Let Em " and Fm the sets of interior edges 
and faces respectively, and accordingly Eqm and Fqm the sets of boundary edges 
and faces. 

The Hopf tensor trace associated with the graph F is the function ttr/^ : 



Q^x.Q^fp Q^y, Q<)g,^ttTr{{xeh{fph{ye};{gg}) 

eeEj^,o peFM eeEoM qeFgM 



(147) 



which is defined via diagrams and the evaluation rules (139), (1401, (141), (1431 
and (fl44|. 



Note that this generalized Hopf tensor trace reduces to the Hopf tensor trace 
for a surface without boundary (as in Definition ^ in the natural way by re- 
stricting any virtual loop function fp to be cocommutative and setting: 



ttrr({xe}; {fp}) :- ttrr({a;e}; {fp}; 0; 0). 



(148) 
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4-2. Quantum states. So far we have merely defined a particular, fully contracted 
tensor network (which is a complex number) with no reference to a quantum 
many-body state whatsoever. We now take the next step and use the tensor 
trace above to generate actual quantum states in a remarkably straightforward 
fashion: 

Definition 4 (Hopf tensor network state). Let Xe,ye G H and fp,gq G X 
as in Definitions^ Let F the graph embedded in the surface M . 

l.IfdM^% then 

IV'r({a;e};{/p};{ye};{5j)) 

:= E E ttrr(K'};{/,};{y^'};{5j) (g) |x^) (g) |y^) . (149) 

2. IfdM = % then 

\i'r{{xeh Up})) ■■= E tt^r({<}; ifp}) (8) l^e) ■ (150) 

eeEM 

In both cases we call the resulting state a Hopf tensor network state on the 
graph F. 

Given an arbitrary region R within the surface M it is straightforward to 
partition such a quantum state into interior and exterior parts accordingly. In 
the language of Hopf tensor network states this simply amounts to cutting virtual 
loops. Without loss of generality we will show this for a graph embedded in a 
surface without boundaries. 

Before proceeding we need to clarify though what a region i? of a surface 
discretized by the graph F is supposed to mean. We will only consider regions 
whose boundary does not cross any edge of F. Equivalently, such a boundary 
may only run through vertices and faces or along edges of F. Any part of the 
boundary which coincides with a graph edge will be called smooth and otherwise 
rough. This coincides with the way a graph may be embedded in a surface with 
boundary (see Figured. It is clear that a valid region R is determined by a 
subgraph Fn up to derormations of rough boundaries, hence we will define a 
region by its associated subgraph. This means that a partition of M into R and 
its complement R naturally divides E into the disjoint sets En and E^. At the 
same time it divides F into the disjoint sets Fr, F^ and Fgn or, put differently, 
boundary faces Egj^ are created which do not belong to either subgraph Fji or 
Fj^. In contrast, there are no boundary edges. Hence the (discretized) boundary 
of the region R coincides with Fgj^ and has a length of |Faij|. By abuse of 
notation we will call this length simply \dR\. 

Proposition 2 (Partitions). Let \ipr{{xe}',{fp})) be a Hopf tensor network 
state on the graph F and R an arbitrary region which is associated to the subgraph 
^R = (Vr, Efj, Fii) C F. Let gq £ X and define 

IV'flXiffJ)) := l^r«(K}fl;{/p}fi;0;{gJ)) (151) 
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where {xe}R — \x^ \ e £ Eji} and {fp}R = {fp \ p G Fr} are the natural 
restrictions to Fr. Then 



l^r({xe};{/p})>= Y. \^M{f',}))^\^R{{0)) 



(152) 



qeFgn 



Proof. Obvious from the complete diagrams and evaluation rule ( 143 ) which 
needs to be applied backwards. D 



Building on this general framework of Hopf tensor network states we can solve 
the generalized quantum double model now. Namely, we identify a particular 
Hopf tensor network state as a ground state of the model. As such this state is 
directly seen to be topologically ordered. Note that we will only make use of the 
structure maps of H as well as the Haar integral and its dual to describe the 
state. Again we assume a surface without boundaries. 

Theorem 3 (Ground state of the generalized quantum double model). 

Let h £ H and G X the respective Haar integrals. The state 



\ipr) ■■= \ipr{h,...,h; (/),...,(!))) 
is a ground state of the D(H)-model. 



(153) 



Proof. Since the Hamiltonian of the D(iJ)-model is a sum of local, commuting 
terms by Theorem [2] it is enough to show that each operator A{s) and B{p) 
leaves the state jf/;/^) invariant individually. In order to do so we may partition 
l^r) into an interior part corresponding to the support of such an operator and 
an exterior part. Both parts are glued via comultiplication in X. It will then 
suffice to prove that this interior part remains unchanged by either A(s) or B{p) 
respectively. 

Hence consider a face p d F with a boundary consisting of r edges. The 
interior part of \ipr) is given by 



J? 



<;^." 



w/i,...,/.))=E'^<nAn>'^ '^ 



i hi, 



(154) 



h'' 



h' 



^ J2 ^iK' ■ ■ ■ K) n f3{s{h';)) \h\) • • • 1^;) . (155) 

with the Haar integrals hi = h £ H and arbitrary /^ G X. It is invariant under 
the action of B{p) as can be seen from 

B{p)\^p{f^,...Jr)) 

= y: m'^ ■ ■ ■ h['^) Hhi'^ ■ ■ ■ /.«) n fAsif^^)) \h^'^) ® • • • ® \h?^) 
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r 

= J2't>'iK'---hr)l[fAS{h'[))\h[)^---^\K) 

= l^p(/l,...,/.)). 

Next consider a vertex s ^ V with r attached edges. In this case the interior 
part of IVt) reads 



\Mh 






K 



h'l 






' hn 



-h[ (156) 



\K). 



(157) 



Note that the orien tation of the graph edges can always be reduced to the above 
setting using (116) for kets and (140) for diagrams. Now we obtain 



A{s)\Mh,---Jr)) 

= E n /^ i^C^'D ^^"+1) E /i (^e^i ) ^2") fr (siK) K) 

(hi)j=2 (h) 

= E n /. (^e^") <+i) E Al^e^i') '^^'''^2") h{s{K) h'i') fr{s{h('^K) K 

{hi)j=3 (h) 

(hi) j=3 (h) 

(hi)j=l (h) 

= \A{fl,---Jr)) 



where we repeatedly used Lemma |4J This concludes the proof. D 
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In fact, the proof of Theorem p^ implies that jf/'r) is also invariant under each 
local action of D(_ff) at any site (s,p), not just under the operators constituting 
the Hamiltonian. More precisely, one has 

Bfis,p) Aa{s,p) IVt) = e{a) /(l^ ) |Vr) (158) 

for any f (S) a G D{H). This can be easily seen from the local D(_ff)-niodule 
structure, the fact that A{s) = Ah{s,p) and -B(p) = B^{s,p) leave It/jp) strictly 
invariant and the properties of the Haar integrals. In other words, the quantum 
state {ipr) is nothing hut a trivial representation of the quantum double T){H). 
Comparing with the comment after Theorem [2] one realizes that {ipr) should be 
viewed as the spatially distributed version ofthe integral 4> ® h € D{H). For 
these reasons one may call {ipr) the vacuum of the model and as such it has 
trivial topological charge everywhere. 

Let us emphasize again that one really needs just a single datum, i.e. the 
finite-dimensional Hopf C*-algebra H, to produce this topological ground state 
since the Haar integral h (and the Haar functional </>) is uniquely defined. In 
particular, the construction is fully basis- independent. 

Furthermore the construction is symmetricT] in the algebras H and X with 
their respective integrals h and 4>, hence it has a natural dual notion. In fact, 
this foreshadows electric-magnetic duality as shown in [T^ . 

While we excluded surface boundaries explicitly for Theorem [3] the following 
example shows what one can learn from the presence of boundaries imposed by 
the underlying surface. 

Example 1. Consider the graph 7^ underlying the diagram shown in Figure [5] and 
let H = CZ2. Then the state 

IV'o):=IV'r({M;{0};{M;{0})) (159) 

is a codeword of the surface code defined in [9j . It encodes the logical state | -I-) . 
The other codeword jV'i) can be obtained by acting on this Hopf tensor network 
state with the appropriate string operator connecting the two rough boundaries. 
This generalizes to any finite-dimensional Hopf C* -algebra by using the appro- 
priate ribbon operators. 

Remark 1. From the proof of Theorem [3] it is clear that any Hopf tensor network 
state of the form 

\^r{K...,h-h,...J\F\)) (160) 

is invariant under all vertex operators A{s). From the perspective of lattice gauge 
theory this means that deforming the ground state \i'r) of the D(i?)-model by 
changing the functions {(/)} 1— >■ {/i} only will never break the gauge symmetry. 
These deformations in the state might therefore well correspond to a local per- 
turbation of the Hamiltonian and thus preserve topological order provided the 
strength of the perturbation is limited to a certain finite threshold [Hl[7] . Some 
initial work towards this direction has been conducted in [14_ which is concerned 
with tensor network deformations of the toric code. 



^ Strictly speaking, symmmetry holds up to a flip in the comultiplication, which is precisely 
the difference between H* and X. 
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4-.3. Hierarchy. Apart from ground states of the D(i/)-niodels the franiework of 
Hopf tensor network states based on an arbitrary Hopf C* -algebra H comprises 
more intriguing examples of quantum many-body states. These can be obtained 
from certain other choices of elements in H and X which are motivated both by 
the algebraic structure itself as well as by ideas of charge condensation 



Definition 5 (Hierarchy). Let A <Z H and B C X Hopf subalgebras with Haar 
integrals hA ^ A and (pB & B respectively. Set 



1^ 



A.Bk 



\Mh, 



l-AifPB, 



■,M)^ 



(161) 



Obviously, with the choice A = H and B = X we recover the state \ipr) 
which is topologically ordered as a ground state of the D(iJ)-model. On the 
other hand, if B = {Ix} is the trivial Hopf subalgebra of X then the resulting 
Hopf tensor network state is a product state for any Hopf subalgebra A C H. 
This can be seen from (23 1. Indeed, suppose a face p has r edges in its boundary 



and fp — Ix- Then iterating ( 23 1 yields 



E 



(1) 



'/^' 



(r) 



LX 



(162) 



and consequently one has 



ttrr({xe};l 



X: 



ix 



n 

eeEM 



(163) 



up to a possible application of the antipode S on each graph edge depending 
on its orientation. Since Ix — ^"'"(Ic) is the analogue of (18) for the Hopf C*- 
algebra X we deduce 



J2{e^ac),S{x'J}{e^{lc),x':) 



(^e) 






siSix'J) e{x':) = <Xe) 



and therefore 



ttrr({a;e};l 



X, 



Ax) == Y\. ^^^"^ 



(164) 
(165) 



eeEn 



which no longer depends on the orientation of F. Now by virtue of ( 150 1 the 
quantum state 



l^rihj 



,hA;l: 



lx)) = 



E 



e{h- 



A.e, 



\h!A..) = 



\hA) (166) 



-eEj, 



eG-Eji 



is seen to factor into a simple product state. 

In between these two extremes a hierarchy of quantum states unfolds which 
are indexed by different choices of A and B. However, depending on the Hopf 
algebra H in question the interior of this hierarchy may collapse partially. This 
means that different pairs {A, B) of Hopf subalgebras may actually define iden- 
tical quantum states. Unfortunately we do not know how to characterize the 
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surviving equivalence classes of states in closed form without additional assump- 
tions on H . 

However, if iJ = CG we have a clear picture of the above hierarchy. It turns 
out that the classes of Hopf tensor network states emerging from the partial 
collapse are isomorphic to ground states of certain quantum double models based 
on groups smaller than G. Indeed, the relevant Hopf subalgebras in this case 
exactly read A — CK and B = £p/'^ where K <Z G \s & subgroup and N <\ 



G a normal subgroup (see Section 2.4.3). For simplicity we abbreviate such a 



pair of Hopf subalgebras by {K, N). Then it is not difficult to see that both 
[K, N) and [K, KON) yield identical Hopf tensor network states. Furthermore, 
if |fci, . . . , k\E\) with kf, & K for each edge e £ E is a basis state then one has 



{ki,...,kiE\\ipr''^^^) = {kik, ■ ■ ■ ,kiE\l\E\\i'r"'^^'^) 



for arbitrary elements l^ E K O N, so all amplitudes are actually constant on 
the cosets ke{K D N). This means that one may apply the canonical projection 
tt: X — > K/{K n N) at each edge and regard (pKnN as the Haar integral of 
£^K/{KnN) ^ Hence the resulting state coincides with the ground state jV'r) of the 
quantum double based on the group algebra of the group K/{Kr\N) ~ KN/N. 
We conclude that the hierarchy of Hopf tensor network states arranges quan- 
tum double models based on different groups in one coherent picture. These 
groups are precisely isomorphic to KN/N so the equivalence classes of states in 
the hierarchy are the trivial representations of the quantum doubles Y){€.{KN /N)') . 
These comprise all possibilities for the residual symmetry algebra after the full 
D(CG)-symmetry has been partially broken down to a smaller symmetry al- 
gebra (4ji5 . For that reason the quantum states in our hierarchy realize the 
condensation of topological charges within the same underlying Hilbert space. 



4.4- PEPS. At this stage it is important to make contact with one of the usual 
formulations of tensor network states. In the PEPS approach |43j a quantum 
state is represented by choosing some fixed basis and encoding the wavefunction 
amplitude for each basis element in a fully contracted tensor network. Contrast- 
ingly, our Hopf tensor network states are defined without reference to any basis. 
It is rather the choice of particular, often canonical, elements {xe\ C H and 
{/p} C X which determines the properties of the quantum state. If the need 
arises one can still obtain an explicit wavefunction expansion in a straightfor- 
ward manner. As the next example shows, those distinguished elements naturally 
generate all relevant sums and amplitudes. 



Example 2. Let H — CP . We would like to explicitly construct the Hopf tensor 
network state jf/'r) for the graph E shown in Figure p^ and obtain an expansion 
in terms of the canonical basis {5g \ g E G} of C*^. We assume the underlying 
surface to be homeomorphic to S"^ as indicated in Figure^ Using the Haar inte- 
gral hi = de of the dual group algebra C*^ and its corresponding Haar functional 
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Fig. 6. Hopf tensor network state \il>r) on a small graph 7^ embedded in S^. The outer circle 
is identified with the north pole of S^. 



G H we get by Theorem [3] 



{h. 



Y,{<l^(^^,^ih'IK))\K)®\h',) 



(h. 



= ^ (l){Su-lSy-l)\Su) <»\Sy) 



u,veG 



k E \^9) ® \Sg 



\G 



geG 



Interestingly, while the Haar integral 6^ merely represents a single element of 
the basis, it is actually the act of comultiplication which produces the sum over 
the entire basis from 6e- □ 

As outlined in the introduction PEPS are defined in terms of local tensors 
with open virtual indices, for a particular example see Q. Contrastingly, our 
diagrammatic notation as given by ( 143 1 uses the comultiplication of the Hopf 



algebra to disconnect and separate local objects on the virtual level, hence it 
appears there are no virtual indices at all in our formalism. 

This is not entirely true. If the functions {/p} for evaluating virtual loops 
belong to a particular clasa^ then indeed one may regard the evaluation of a 
virtual loop as tracing over a product of certain matrices, once a particular basis 
of H has been chosen. Consequently, each such matrix naturally constitutes part 
of a local tensor with open virtual indices and evaluating the loop corresponds 
to a cyclic contraction of those indices. Note that this applies to all Hopf tensor 
network states of our hierarchy (but is not limited to these). For example, for 
the state \ipr) at the top of the hierarchy one has the local projector 



P = 



hiri: E (4^'l,«')..l^')(a,/3,7,'5| 



(167) 



{h) a.,l3,-i,5eB 



* This is the so-called character ring Rc{H) 
characters of H. 



X]"=i Cxi where the Xi a^re the irreducible 
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where (L^)q^ denotes a matrix element of the action ( 102 1 with respect to some 
basis ;B[j Note that in this setting the local matrices encode both information 
about the spin state on the graph edge as well as about the function used to 
evaluate the virtual loop. 

Furthermore, if a Hopf tensor trace admits a representation in terms of local 
tensors with open indices properties like injectivity may be studied. More pre- 
cisely, a PEPS is called injective [3S] if there is a partition of the underlying 
graph in disjoint regions Ri such that for every region Ri the linear map from 
open virtual indices at the boundary to physical indices in the interior is injec- 
tive. Whether or not a particular tensor network representation of a quantum 
state is injective has important consequences for e.g. the existence of a parent 
Hamiltonian that has the given PEPS as its unique ground state [32]. One can 
show that all Hopf tensor network states {ipp ) in the hierarchy are not injective 
for any finite-dimensional Hopf C*-algebra H. However, (at least) the state {i/jr) 
at the top of the hierarchy obeys a relaxed version of injectivity which one might 
call H -injectivity. This is a certain generalization of the G-injectivity condition 
defined in [iTj . 



5. Calculating the topological entanglement entropy 

Topological order is commonly associated with non-local order parameters. It is 
believed that among these resides the topological entanglement entropy 7 fI7\ 
I32j which is a universal additive correction to the area law for the entanglement 
entropy of a bipartition of the system into a region R and its complement. Given 
a ground state of a topologically ordered system, the entanglement entropy Sji 
for a region R is argued to scale as 

SR = a-\dR\-j (168) 

in the limit of infinitely large regions. Here a is a parameter which encodes 
non- universal behaviour on short length scales. In fact, the topological entangle- 
ment entropy 7 also contains partial information about the types of quasiparticle 
excitations that may occur in the low-energy sector of the system. 

For these reasons we would like to show in this section how the non-local 
order parameter 7 can be understood naturally in the context of Hopf tensor 
network states. In particular, we are going to show that different classes of finite- 
dimensional Hopf C* -algebras yield fundamentally different mechanisms for the 
emergence of a non-vanishing topological entanglement entropy. 

To this end, we will exploit the generic decomposition of Hopf tensor network 
states into interior and exterior parts as stated in Proposition |2] We will then 
compute the block entropy directly from certain properties of the reduced density 
operator. In order to render its simple structure evident we will make use of 
isometrics and completely clear out the interior of both the region R and its 
complement R. Effectively, we concentrate all topological information contained 
in the ground state into the (inner) boundaries of the system. Furthermore, the 
distillation process will be crafted such that all intermediate quantum states can 
be kept track of conveniently via Hopf tensor traces. 



^ Furthermore one may express the (physical) ket in the same basis, too, and appreciate the 
similarities and differences as compared to the trivial case (|10|. 
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It should be noted that our scheme of applying isometries implements entan- 
glement renormalisation |^H5] both for the D(_ff)-models as well as for states 
of the hierarchy. In fact, it directly extends the work in [1 which is concerned 
with states at the top of the hierarchy for H = CG. At the same time our scheme 
provides a complementary view on entanglement renormalisation for string-net 
models [30] , namely from the perspective of a local symmetry algebra. 



5.1. Isometries. We begin by developing the distillation process. For that we 
will introduce a hierarchy of "little" isometries that insert minimal faces and 
vertices into a given graph. The little isometries at the top end of the hierarchy 
are taylored to the generalized quantum double model, i.e. they preserve all 
excitations of that model. Subsequently we will define a pair of local unitary maps 
that allow for reconnecting edges. Taken together, both the little isometries and 
the unitary maps will yield a hierarchy of isometries that add or remove arbitrary 
edges. Their effect on Hopf tensor network states will be analyzed in the next 
section. 

As far as notation is concerned we will always denote the original graph by 
Fi and the modified one by /Ij. Whenever it is appropriate to talk about actual 
quantum double models, Hi will denote the Hamiltonian of the D(_ff)-model on 
the graph Fi. 

Proposition 3 (Little isometries). Let A C H and B C X Hopf subalgebras 
with the respective Haar integrals Ha G A and (pB & B. Additionally, let 

(169) 

(170) 
(171) 



Then the linear maps ip 



^A,B 


■.^Y.^B{h':^)^'A^ 




[hA) 


^A,B 


^A.B 

II \ II ' 

\\-^a,b\\ 


Aa 


_ Ha _ 1 




\\hA\\ ^4>{hA) 


^A,B 


and iy defined via 




^9 + 1 



(172) 



and 




Xg Xg-I-X 




(173) 



Xq Xq^;l 
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are isometries. 



Furthermore, for ip :~ ip' and iy '-^ iy one has 



ipHi = 'H2iFi 
iyHi ~ T-L2W- 



(174) 
(175) 



Proof. For the first claim, it is enough to show the invariance of the inner product 
or equivalently that ||1/i.b|| = W^aW = 1- Indeed, from 

\\hA\\^ = {hA, hA) = (f){h*^hA) = 0(/ii) = H^a)- 

this is easily seen to be true. 

As for the second claim we note that ip inserts a loop which is automatically 
stabilized by the corresponding face operator in 'H2- It is also easily seen that 
1h,h' = Iff and hence it does not affect the magnetic flux through the original 
face as measured by either "Hi or 'H2- Since 1^ is trivially invariant under the 
adjoint action we see that the vertex operator at the leftmost vertex is not 
affected by the additional face, either. By similar arguments one proves that iy 
intertwines the Hamiltonians, too. D 

For the sake of completeness we also give the maps (ip' Y and (jy)^: 




Xq + l 




{4-^)' 



> (U,B,fc) 



(^^)' 



>{AA.b) 




Xq Xg_|_l 



Xq Xq^\ 



(176) 



(177) 



Now we define the unitary maps that allow for reconnecting edges of the 
underlying graph. 

Proposition 4 (Unitaries). The linear maps Up and Uy defined via 




(178) 



Xq + l 
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and 



6<i'a;i 




Uv. 



E 



Xg 3^g_|_l 

are unitary and their inverses read 



b'^'a;. 



^g+i 



(179) 




Xq + l 



where c — S{x' ■ ■ ■ x[)b and 



S{b^'''>)xi 



{,(9+1) 



X, + l 




respectively. 

Furthermore one has 



S(6<i')a;, 



UpTLi — 1-L2Ufi 
UvHi =H2Uv 



Xq + l 



(180) 



(181) 



(182) 
(183) 



Proof. It is easy to check that Up (Uy) as defined above is both a left and a 
right inverse of the map Up (Uy)- Indeed, leaving out the part x^+i • • • (X) x^ 
(which is unaffected by Up) one has: 



' x'' <Si x' ■ ■ ■ x'lb 



ulUpixi (g) • • • (g) x, (g) 6) = C/U ^ x" 

= ^ < (g • • • (g x'^' (g S{x'^ ■ ■ ■ x'l) x; • • • x[b 

(Xi) 

= E<®- 

= ^a;'/(g- 



'x'^'^Six';)---S{x'^)x'g---x[b 
Xg g) e{x[) ■ ■ ■ e(x^) b 



{X^) 

Xl (g • • • g) Xg g) 6. 
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Similarly, one shows UpUp ~ id. Now 

(b) ((,(9 + 1)) 

^5((5(?+i))(i)) 6(9)3.^ ^(5('?+i))(«+i) 
= J2 S{b'^^'^^) b^^hi (g) • • • ® 5(6(«+i)) b^'^^xg (E) &(29+i) 

= ^ S{b") b'xi (S)X2<»---<»Xq(» b'" 

— X\® ■ ■ ■ ®Xq®b. 

and by the same token one proves UyUy — id. 

For the first claim it remains to show that for all Xi,yi^b,c G H the inner 
product (|119p is invariant under Up [Uv)'- 



{Uf{xi (g) • • • «) a;, (g) 6), Upivi g) • • • g) % «) c)) 

= f ^ x'/ (g • • • g) a;^' (g x^ • • • a;'i6, Y^Vi ^ ' ' ' ^ Vq ® Vg' ' ' Vi'^) 

= J2 (.x'^---x[b,y'^---y[c)f[ix';,y';) 
(a;i)(yi) J-i 

g-l 

(a;i)(yi) i=i 

9-1 
= (x„2;,) ^ </'(&*(a;D'---K-i)X-i---yl^)n'^((^P"y") 

= (xi g) • • • g) Xq g) 6, 2/1 (g • • • (g y, (g) c) . 



Note that we used property ( 38 ) of the Haar integral cj) in the fifth line. Further- 
more we have 

{Uvixi g) • • • (g Xq g) 5), [/^(yi (g • • • (g yq (g) c)) 

= ( E ^'^^^1 «) • • • g) b^'^'>Xg (g 5(«+i), ^ c(i)yi g) • • • (g c(«)% (g c^-J+i) I 
^ (ft) (c) / 
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(fc)(c) i=i 

(bKc) 3 = 1 

(6*c) i = l 

9-1 

= K,y,) 5] cj^{ib*c)^'^^)Y[<j>{x*{b*c)(^k,) 

(b'c) j = l 

q 

= ib,c)Y[ixj,yj) 

j=i 

= (xi (g) • • • (g) a;^ (X) &, yi (X) • • • (g) 2/g (X) c). 

For the second claim, one needs to show that Up (Uv) commutes appropri- 
ately with the terms A{s) and B{p) in the Hamiltonians Hi- This can be done 
directly for those faces which are complete in (178). Stating this more precisely 
we have: 

UfB{p,)^B{p',)Uf ie{i,2} 
where faces are labelled as follows: 





Indeed, consider 

UfB(p2){xi (E)---(E)Xq(E)b) 



Upixi 



X, (X ^ (f,{b') b" 

(b) 

' Xg (E)Xg- ■ ■ x[b'(l){b") 



(^.)(fc) 



= Y, ^{S{b") 5(xf )) • • • S{xf^)xf^ ■ ■ ■ x^^^) x['^ ®---®xf^® xW • • • x^pb' 
(^.)(fc) 

. Y, 0(^(4'^ • • • ^f^b") xf • • • xf^) x'^l^®---® xW ® xW • • • x^pb' 
(^.)(fc) 

= (id (X • • • (X id X)S') ^ Y H^i^'q ■ ■ ■ 2^'i&)'< • • • x'l) 

{Xi}{b) {S{x'^-x[b}) 

x'i"X)---(Xx;"X)S'(x;---x'i6)" 



A hierarchy of topological tensor network states 



45 



= B(p2) Uf{xi (E)---(S)Xq(E)b). 

Note that in the fourth hne we repeatedly used the definition of the antipode. 
Hence both B{p2) and i?(p2) project h onto the value (j){h) 1h- In the following 
we may therefore assume that b is replaced by 0(6) 1h- We now analyze the 
other complete face: 

UfB{pi){xi ® ■ ■ ■ ®Xq® (p{b) 1h <8) Xq+i (X) • • • (g) Xr) 

= <j){b) Up i ^ (jjix'r ■ ■ ■ x\)x'{ ®---®x'q®lH® X^'+i (g) • • • ® x" j 



x'q ®x'q--- x'-l (X) x"4.i (g) • • • «) x; 



{^^) 



1 yy ^g+l ' 

^ Xq (S {x'q ■ ■ ■ X[y 



= Hb)Y. E '^(4---^;+iK---^'i)')^i®- 

(a;i) (a;;j---a;'i) 

= 0(6) S(pi) f ^ x" (g) • • • (Xl X^' «) 2;^ • • • x'l (g) Xq+i «) • • • (g Xr j 
= B{pi) Uf{xi ® ■ ■ ■ ® Xq® (j){b) \h ® Xq+i ® ■ ■ ■ ® Xr)- 

For those surrounding faces which are affected by Up one may easily prove 

[Uf,{tI),] = q ze{i,...,g} 

for any f <E X. Here (Tjf )i deno tes t he action on the i-th tensor factor. 

Those incomplete vertices in ( 178 ) which are affected by Up can be dealt with 



if the following is true for all a G i/, i € {1, . . . , g — 1} and b ex 1h'- 

UpiL't)! = fe(i'i')i ® {L't'h)uF 



(a) 



C/F,^(i^'),®(L«"),+ i 



(a) 



= 



C/f ^(<),0ad(a% = ^(if ),®(i;')o C/^. 

Here denotes the bubble edge that is reconnected by Up- While it is easy to 
show the first two equations, the remaining one needs some more attention: 



Up i ^{Ll)q (g> ad(a')o j (xi (g) • • • (g x, g) 0(6) Ir) 

= Up(^xi'^---'^ Xq_i (g a"xq g) 0(6) ad(a')(lH) 
= 0(6) Up{xi g) • • • g) Xg_i g) axg (g Ifl-) 
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• (E) Xg_i (g) a"x'^ (g) a'x'g ...x[ 






E(^+")<? ® iL+)o jUpixi g> ■ ■ ■ (g) Xg (g> (t){b) Ih) 

(a) ^ 



This completes t he p roof of (182 1. 

The proof of (|183) is analogous and left as an exercise for the reader. D 



Corollary 1 (Isometries). The linear maps Ip' :— Up ° ip , 




-E 




(184) 



and ly :— Uy ° «y, 




A^"'^,, 



Xq+l 



are isometries. Furthermore, for Ip 



Ip' and ly '■= ly one has 



(185) 



lyUi - n2ly. 



(186) 
(187) 



Put differently, the collection of these maps allows to move between the D(iJ)- 
models defined on two arbitrary graphs which are embedded in the same surface. 
Furthermore, all maps are local, i.e. they cannot change nonlocal topological 
quantum numbers. 



5.2. Transforming hierarchy states. In the preceding section we have learned 
that the set {If, Iv} of isometries maps the ground state subspaces of generalized 
quantum double models on related graphs onto each other. We would like to 
prove that this not only holds for the entire subspace but also for individual 
ground states. Namely, the isometries {Ip, ly} precisely identify the Hopf tensor 
network states \4'ri) as given by Theorem [3] for the different graphs Fi. 

Keeping our focus on states we will furthermore show (for a special case) that 
the hierarchy of isometries {Ip ' ,Iy} has the very same effect on the hierarchy of 
states {-tpp^ ) as given by Definition 5 In other words, identifying these states in 
an isometric fashion is as simple as adding or removing edges from the underlying 
graph and adjusting the tensor trace canonically. 
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Once we act with the isometries on our hierarchy of Hopf tensor network 
states it will become clear that the quantum states defined in Section |4.2| are 
not properly normalized. We can (partially) remedy this situation by choosing 
integrals of unit norm instead of the Haar integrals Ha and (J)b ■ In other words, 
we need to include a factor of ||ft.A||~^ for each edge and a factor of ||0b||~-^ for 
each (complete) face of a Hopf tensor network state. It is understood that from 
now on all Hopf tensor network states will be normalized in this fashion without 
change of notation unless otherwise stated. 

We now fix some notation in the following definition so that we can easily 
refer to pieces of Hopf tensor network states later on. 



Definition 6. Let A C H, B C X Hopf subalgebras, h^ 
respective Haar integrals and let fi, . . . , fr '^ B. Set 



e A, 



e B the 



fa 'i^-" 



(hA.i) 



/a+i 





(188) 



I V-pup' ({/<})):= 



\hA\V+'\\M\ 



E 



{hA,i)(hA) 





(189) 
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te"((/.})) := P^ E 

('lA,i) 





^A,q + 1 



(190) 



and 



" " {hA.^){hA) 



h\ 



h\ 




/g+l 



/« ''A,(j+1 



(191) 



/j'. 



''■:4.,+i 



As before we simply write |V'r(/i, ■ • ■ ,/r)) whenever A = H and B = X. 
Note that in this case one deals with a fragment of a quantum double model 
ground state! 



Proposition 5. 



I F\lpp{fl, ■■■,!,■)) = \ll'pUp'{fl,---Jr)) 



/i?IV'f'^(/i,-..,/.)) = ICuf'(/i. •••./.)) 



(192) 
(193) 



Proof. First let us prove the relation involving the face isometry. It is enough to 
consider a face which is bounded by two edges. In this case we first transform 
the right hand side of the equation in order to eliminate one of the virtual loops: 

2 

|V'pup'(/i,/2)) = \H\'^ Y. <l>{hrh'")Hh2S{h"))l[f,{S{h'^)) \h[) ® \h') ® \h'^) 

{h.){h) 3 = 1 

(In) 
= \H\iY.^{K'h['^) h{S{hf^)) h{S{K)) \h^?) ® \S{h^i^)) ® \K) 

ih^) 
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where in the second hne we employed Lemma [5] 
In diagrammatic notation the above amounts to 



|V'pUp'(/l,/2)) = |i?|'E^'^<l 




(hi) 



which proves the claim about the face isometry. 

As for the vertex isometries we can in fact assume A ^ H without loss of 
generality. By the same token we will never make use of fi E B in the following 
argument, hence it suffices to assume B = X, too. We begin with the state on 
the original graph: 

r 

\Mfi,---jr)} = \\h\rJ2Y[f,{s{h';)h'f^,)\h[}^---®\K). 
{h.)j=i 



Applying ly yields: 

\\h\\-+'lv\Mh,---Jr)) 

r q 

(h)(lH)j = l fc=l 

= Y, fr{S{K)h'[')f{f,{S{h';)h';'^,)<^\h(^^h',)®\h('>+^^)®--- 
(hi)(h) j=l k=l 

where in the last line we have abbreviated the remaining tensor factor 



n fAs{h/;)h/;',^ (g) \K) 

j=q+l k=q+l 

by the trailing ellipsis. In order to simplify notation we will continue to do so in 
the following. Hence: 

||/l|r+l/y|V.(/l,...,/.)) 

= ^ fr{s{K)hr)h{s{h';)h'^')f[f,{s{h'^)h'^'_,,) 

{h,)ih) J=2 

l/i^^^/i'i) ®---® \h^'^^h'^) ® |/t(''+^)) ®■■■ 

= ^ U{s{K) s(h(^y) K) h{s{h'i) /.(2)/,-) f,^s{K) h'i') n f,{s{h';) /.- 1) 

(lu)(h) i=3 
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{h,)ih) j=3 

\h[) ^ \h'2) ^ 1/1^^^/13) (»•••(» Ih'^'^h'g) ® \h'-^+^^) ® • • • 

= Y^ fr{s{K)s{h('->)hr)f,{s{h';)h^'^h';'^,)l[f,{s{h';)h';'^,) 

{h.){h) j=l 

\h[) ^ ■ ■ ■ ^ \h'g) ^ \h^^^) ® ■ ■ ■ 

= E Us{h';)h'''h'^'^,)fr{s{K)s{hnhr) l[fMh7)hr+i) 

(h,)(h) j^q,r 

\K) «) ■ • • ® \h'q) (g) \h') (g) |/l^+i) (g) • • • ® \h'^) 
= ||/l|r+l|V^.U.'(/l,...,/r)). 

Note that from the third hne on we repeatedly applied Lemma |4j D 

Unfortunately, we do not know how arbitrary hierarchy states {ipp' ) are 
affected by the associated face isometries Ip' for a general finite-dimensional 
Hopf C*-algebra. However, for H = CG we can state the following and leave the 
proof to the reader. 

Lemma 3. Let K G G a subgroup and N <i G a normal subgroup. Furthermore 
let fi e C<^/^. Then 

/^■^ l<^~(/i, . . . , /.)) = ^/^ |^^^^;(/,, . . . , /,)) . (194) 

The above lemma shows that unless one draws the quantum state from the 
top of the hierarchy one needs to include an additional factor of ^y\KN\ / \G\ 
for each face to insure proper relative normalization of the Hopf tensor network 
states. In order to fix the normalization absolutely we may calculate the norm 
of a single (preferably small) Hopf tensor network state on a given surface. 
Since all other states on the same surface can be reached from this initial one 
via isometries their norm will be determined automatically. It turns out that the 
absolute normalization factor is entirely a property of the surface the Hopf tensor 
network state is embedded in. For example, absorbing the additional factor per 
face into the definition of the Hopf tensor network states one has for arbitrary 
graphs r on 5^ 

m'^'''}\\^^/\KN/N\ (195) 

while on T^ one has 



Here Ccig) denotes the centralizer of the element g in the group G. In order 
to keep the following discussion as general as possible we will stick to relative 
normalization unless otherwise noted. 
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5.3. Entanglement entropy. Having defined the isometries Ip and ly we finally 
embark on an exact calculation of tlie entanglement entropy for a simply con- 
nected region R on S^ . We focus on the Hopf tensor network state IV'r)- 

In the following we show how the inner product between two face pieces of the 
Hopf tensor network state \4'r) depends on the boundary degrees of freedom. 

Proposition 6. Let fi,gi G X and r — \dR\. Then 

r 

{i;p{h, . . . , fr)\M9u- ■ ■ ,5r)) = \Hr E '^('^1 • ■ ■ ^^) U^ajflKf^'D- (197) 

(hi) J=i 



Proof. Let Aj = hi = h the Haar integral of H. Then from ( 155 1 we get 



(Mfu- ■ ■ . fr)\MgU • • • ,5r)) = J2 -^e^"' ■ ■ • '^l") '^(^"' ■ ■ • ^i") 



= Y^ cj^{h'i'---K')^{K'---x'i') 

(h.)(A,) 
r 



where we used (119), the properties of h and the involution * as well as (51 1 for 
the second line. Employing Lemma |6] we can simplify this to 

(Vp(/i, . . . , /.)|Vp(5i, ■ ■ • ,5.)) = \H\-' Y. ^(^"1 ■ ■ ■ K') ^i^r ■ ■ ■ ^2 S{h[)) 

r 
(hi) 

Finally, the claim follows from (jy^ — (p. D 

Analogously, we are interested in the inner product between vertex pieces of 
the same Hopf tensor network state. 

Proposition 7. Let fi,gi G H* and r — \dR\. Then 

r 
{Aifl, ■ ■ .JrMsigi, ■ ■ .,9r)) - l^fE n ni^fSihfl)) g,{hfS{hfl{)). 

(198) 
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Fig. 7. Boundary graphs for a simply eonnected region R C S^. Here black edges constitute 
the corresponding subgraph Fn while grey edges belong to F^. These graphs can be reduced 
further (without touching the boundary) to minimal graphs. 



Proof. Let again A,; = hi = h the Haar integral of H. From (157) one has 



-EEn^(^"^(^^i))5.(5(A^)A^;0'^('^^A^)- 

(h.)(AOi=i 
Additionally, Lemma [6] yields 

Eri<;4(a.A;.)5,(5(A;') A- = li/fEn^.K'^K+i)) 

(Ai)j=l (a.)J = l 

where aj e iJ is arbitrary. Armed with this identity it is easy to verify the claim. 

D 

We now state the main result of this section. 

Theorem 4. Let R d M a simply connected region and let Fa C F the corre- 
sponding subgraph. Then for any a > the state lipr) has the Renyi entangle- 
ment entropies 

Sa{pB.)^\dR\log\H\~log\H\. (199) 

In particular, the topological entanglement entropy reads 



7 = log \H\ 



(200) 



Proof. Without loss of generality we restrict to M ~ S"^ . Acting with the isome- 
trics If and ly on \ijjr) we may wipe out the bulk of both R and R completely 
and reduce F to one of the graphs shown in Figure [7] where we indicated several 
types of boundaries. 

We further assume the first type of boundary in the following. In this case 
we can continue the isometric reduction to the graph Fq given by 



1 R 

2 



and for the resulting Hopf tensor network state on this graph Proposition |2] 
yields the natural splitting 



l^ro) 



E I V'«('/''i,...,<^;))®iV'fl('/'i, •••,'/'")) 
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where the states 

\i,R{h,...Jr)}^\\h\rM-'\Mfl^---Jr)) 
\i^R{fl,...Jr))^\\h\rHsifr,---Jl)) 



are identical to (154) and (156) respectively up to normalization. Note that 
r — \dR\ denotes the length of the boundary. 

We determine the reduced density operator p^ — N'^pR of the region R as 
follows. Let (pi = (fi = (j) the Haar integral of X and hi = Xi = Ki = h the Haar 
integral of H. Ignoring the normalization factor N for a moment one has 



PR = ^^RiHro) (VToI) 

= J2 l^«(<^i' • ■ • ' ^'r)) i^Riv'l, ■■■,Vr)\- {^Ri^'l, • • ■ , ^rMRi^l^l, • ■ • , €)) 

{4>.){v>.)^{\,) fc=i 



5]</.«"...on¥'K5K')) («i,---:4l 



{h,)j=l 

[K) («>) (A.) J = l 

(l){hfs{hfl^) <) \X\) ®---® \K) ® 0(4?) ®---® <^«?) 



where we used Proposition [7] in the third line. 
Furthermore, using Lemma [7] we may derive 



E '^« • • • <") n '^(^^4) (8) '^(4?) = i^r'' E '^(^" • • • ^i' ) ^is%) ?) 

(k.) J = 1 fc=l (b.) .7 = 1 
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With this under our beh we arrive at 

n <i>{hf s{h^},) s{x'^)) \x[, . . . , a;) <S)<t>[s{hi'^s{h^^l,)) ?] 
j=i k=i 

= \H\ J2 HK' ■ ■ ■ A'l") n Hhf^Sih^"^!) si^'D) 

r 

fc=l 
Again some more preparation, namely Lemma IS] allows us to find 

r r r 

(A.) i = l fe=l (a,-) J = l 

and subsequently we have for the reduced density operator: 

(h,) ]=i k=i 

Finally, it can be shown that the reduced density operator of the region R takes 
the simple form 

pn = Iff r-^ E I'^i'^) ® ■ • • ® I'^-i) ^ 1^(^-1 ■ ■ • ^i'^)) 

C*.) (201) 

® (l){Sihf'>) ?) ® • • • ® 0(5(/if_\) ?) ® <^(/i|,i\ • • • /if ^?). 

up to normalization. It is easy to see that 

and hence we can fix the normalization by setting iV = |iJ| 

Furthermore it is not difficult to show that pji is proportional to a projector. 
Indeed, consider 

® (/.(^(Af ^) 1)®---® (/.(^(AfJi) ?) ® ,/)(a|.'Ji • • • A^'^?) 

= \Hr'PR. 
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This means that the spectrum of pj^ is flat. 

Finally one can prove that the rank of the reduced density operator pn is 
given by |i?|'^~ hence each non-zero eigenvalue equals \H\ ^^ . Then the Renyi 
entropies read 

^"(z^^) = 13^ log(tr(pS)) = \dR\ log \H\ - log \H\ . (202) 

independently of a. D 

Remark 2. We would like to stress that any reduced density operator pn obtained 
from the states \'ipK,N) is proportional to a projector as has been noted before 
for the case of group algebras at the top of the hierarchy [2D]. In other words, 
the entanglement spectrum is flat and all Renyi entropies are equal. 



While this proof yields exactly the entanglement entropy one would expect it 
does leave open one question: what is the mechanism of this universal correction 
to the area law in these Hopf tensor network states? 

It turns out that one can learn the answer by analyzing the boundary con- 
figurations {fi,...Jr) C X'' of cither |V'fl(/i, • • • , fr)) or |i/'fl(/i, • • • , fr)} in a 
particular canonical basis which is natural to the problem. Namely, by a gener- 
alized Fourier construction [10,12 one can extend the grouplike elements Q{H) 
of -ff to a linear basis the dual of which we denote B G X. This dual basis is 
essentially unique and therefore serves as our canonical choice. For example, if 
iJ is a group algebra one trivially has B = {5g \ g € G}. 

Then the dependence of, say, the interior states |'!/'k(/ij • • • j fr)) on the bound- 
ary configuration (/i,...,/^) C B"^ is radically different if 7? is a group al- 
gebra, the dual of a group algebra or a non-trivial Hopf C* -algebra like H%. 
Indeed, one may show that in the first case all non-zero interior states are or- 
thogonal and share the same norm. Furthermore an interior state with bound- 
ary configuration {5g-^, . . . ,5g^) vanishes precisely if gi • • • g^ 7^ e. Hence the 
topological constraint is realized by simple group multiplication around the 
boundary. In the case when H equals the dual of a group algebra, any inte- 
rior state |'0/f(/i, . . . , fr)) is different from zero, however, altogether they cease 
to be orthogonal to each other. Instead, they come in bunches of exactly \H\ 
identical states each. This obviously defines an equivalence relation among dif- 
ferent boundary configurations. It is only these equivalence classes which are 



orthogonal and need to be summed over if we want to turn (5.3) into a proper 
Schmidt decomposition. This is the second mechanism how the topological con- 
straint is implemented in a quantum double model ground state. Finally, if H is 
a nontrivial Hopf C* -algebra like H^ one will observe a combination of the two 
mechanisms just described. 



6. Discussion 

First of all, wc have developed a tensor network language, based on the formalism 
of finite-dimensional Hopf C*-algebras, which is both more flexible and more 
natural than conventional descriptions for topological states. We have given rules 
to evaluate tensor network diagrams by means of tensor traces and to construct 
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quantum states on the lattice, as well as to perform operations related to the 
spatial lattice structure, such as cutting and joining along region boundaries. 

This tensor network language has been shown to directly lead to the con- 
struction of well-known topologically ordered states, namely ground states of 
Kitaev's quantum double models based on groups [35]. But in addition, we also 
obtain the extension of Kitaev's construction to Hopf algebras. This is physi- 
cally relevant because the quantum double models based on Hopf algebras is the 
smallest class of models allowing for a non-Abelian electric-magnetic duality in 
topological models |12j . All these states are written, in a basis- independent way, 
in the form of tensor networks involving the intrinsic Hopf C*-algebra structure 
only. 

Relaxing this latter property leads us to a hierarchy of states defined from 
different subalgebras of a given Hopf C* -algebra. We study the classes obtained 
from group algebras, that is, Kitaev's original models, and show that the hier- 
archy arises from the mechanism of condensation of topological charges [HH], 
and so we conjecture that this mechanism can be described in general in our 
language. 

Furthermore the hierarchy states can be regarded as ground states of certain 
frustration- free Hamiltonians beyond the D(if )-model. These are obtained from 
Theorem [2] by replacing the Haar integrals h G H and (j) G X with Ha & A and 
(pB G B respectively where A C H and B C X are Hopf C*-subalgebras. 

Not least, we have established isometric mappings defining entanglement 
renormalization [45ll46] extending the work of [1] for the states at the top of the 
hierarchy. This is a systematic procedure to thin out degrees of freedom keeping 
the topological nature of the states, and hence their nonlocal properties, intact. 
Our computation of the topological entanglement entropy is an application of 
this general scheme. 

From a more philosophical viewpoint, our graphical language brings tensor 
networks closer to algebraic sources; and, in particular, it highlights the role 
of Hopf C* -algebras as natural spaces for transformations and symmetries of 
quantum many-body systems. 

We also give a broader perspective of Kitaev's original quantum double con- 
struction, giving the mathematical background necessary but emphasising the 
physical meaning of the different elements. Hopf C*-algebras are shown to be 
natural and powerful tools in the practical analysis of two-dimensional quantum 
problems; we hope that the applications in this paper constitute a learn-by-doing 
introduction to this structure for physicists. 

This work opens several areas of research we are currently pursuing. First 
of all, it is important to study the precise form in which charge condensation 
appears in the Hopf quantum double models. Next, it is natural to ask whether 
all topological models can be cast in quantum double model form by a suitable 
generalisation of our tensor networks. In (12] we conjecture that the answer is 
yes, and that the relevant algebraic structure is the class of weak Hopf C*- 
algebras; the details of the mathematical construction and the corresponding 
tensor network are work in progress |13j . 
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A. Hopf singlets 

We collect here some lemmas, with their proofs, that are used throughout the 
text. 

Lemma 4. Let h (z H the Haar integral, a,b,c Cz H and f^gCzX. Then one has 
Y,f{S{h")b)g{ch"')ah'^ J2 f{S{h")a"b)g{cS{a')h"')h'. (203) 

ih) (a)(h) 

Proof. It is clear that the statement will follow from 

^ ah' (g, Sih") (g, h"' ^^h' ® S(h") a" ® S{a') h'" . (204) 

[h) (a){h) 

Indeed, we have 

^ ah' (g, Sih") (g, h'" = Y, a''^' ® ^(^") ^("") ^ ^"' 

(h) (a){h) 

= J2 a'h' g, S{h") S{a") a'" g, h'" 

{a){h) 

= Y, a^'^h' g> S{a<^^^h") flW (g, £(a(3)) h'" 

ia){h) 

= Y a(i)/^'®5(a(2);i")a(5)^5(a(4))a(3)/j"' 

{a){h) 
(a) (a'h) 

= Y <«') ^' ® S{h") a'" g> Sia") h'" 

ia)ih) 

= Y h' g)S{h")a" g)S{a')h"'. 



(a)(h) 



D 



Lemma 5. Let a,b Cz H and h Cz H as well as (/) Cz X the respective Haar 
integrals. Then 

Y (^{ah'") (j){b Sih")) h' = \H\~^J2 <^(^"') ^(«")- (205) 

(h) (a) 

Proof. Similarly to the proof of Lemma |4] one can first show that 
^ /i' ® Sih") g, ah'" = Y Sia") h' g, Sih") a' g, h'" 

(h) {a){h) 
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holds for any a ^ H. Then we use the defining property (38) of the dual (Haar) 
integral: 

^ (j){ah'") (t,{b S{h")) h' ^Yl ^i^'") ^(.^ ^(^'") «') -^("'O ^' 
(h) {a){h} 

= Y,^[bS{h"q}{h"'))a']Sia")h' 

{a){h) 

= ^ (f>{bS{lH)a')S{a")h'(f>{h") 

ia){h) 

(a) 
D 

Lemma 6. Let a,b,c,d ^ H and f,g e X . Furthermore let h ^ H and (/) (z X 
the respective Haar integrals. Then 

Y,<P{ah')f{hS{h")c)g{dh"') - \H\-'Y.f{ba"c)g{dS{a')). (206) 

(h) (a) 

Proof. From ( |204| we get 

^ <j,{ah') f{b S{h") c) gidh'") = Y. '^(^') /(^ ^e^") °"c) 9{d S{a') h'") 
CO (")(?») 

(a)(;>) 

= <^(/i)^/(5a"c)5(d5(a')) 

(a) 

which proves the claim. D 

Lemma 7. Let a,b,c,d ^ H and h E H and (p E X the respective Haar integrals. 
Then 

Y4'ih'a)^{bh")^{ch'"d) = |iJr'^0(5(6')a)(/.(c5(6")d). (207) 

(h) (b) 

Proof. We first show that 

Y,h'0 bh" ® h'" = Y S{b') h' (g, h" (g> Sib") h"'. 
W (fc)(/i) 

Indeed, we have 

Y^h'^S) bh" (» h"' = Y <^') ^' ® ^"^" ® ^" 
(ft) (b)(h) 

= Y S{h') b"h' ® b'"h" ® h"' 

(b)(h) 
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ib){h) 
(b)(h) 

(b) (b"h) 

= ^ 5(6') e(5") /i' (» /i" ® 5(6'") h'" 
(b)(h) 

= ^S(b')h' ®h" ®S{b")h"'. 
{b){h) 



Again by property (381 of the dual (Haar) integral we deduce: 



^ (j){h'a) cj){bh") (j){ch"'d) = ^ 0(5(6') h'a) (f>{h") (f>{c 5(6") h"'d) 
(h) {b){h} 

= J2 H^ib') a) <p{cS{b") (f>{h') h"d) 
(b)(h) 

= ^{h)Y,cl,{S{b')a)(^{cS{b")d). 

O 

Lemma 8. Let a,b,c £ H and h (z H as well as 4> £ X the respective Haar 
integrals. Then 

Y, (j>{h"a) <i){bh"'c) h' = \H\-^ Yl ^{bS{a") c) S{a'). (208) 

i.h) (a) 

Proof. Along the same lines as the proof of Lemma [7] one can show that 
^ /i' (g) h"a (g) h'" =Y^' '^("') ® ^" ® ^"' '5'(«")- 

(h) {a)(h) 

As before the rest then follows from the properties of the dual Haar integral. D 
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